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ABSTRACT 


A  theory  (Trof imenkof f ,  et  al ,  1972)  of  a  non¬ 
equilibrium  electron  tunneling  in  normal  metal- 
insulator-metal  junctions  has  been  proposed  to  explain 
the  presence  of  a  small  (^1%)  conductance  reduction 
around  zero  bias.  In  particular,  the  related  struc¬ 
ture  in  a  and  da/dV  have  been  accounted  for  by  the 
blocking  of  otherwise  available  electron  states 
due  to  finite  electron  relaxation  times  in  metal 
electrodes.  Detailed  measurements  on  Mg/Mg,  Mg/Au  and 
Al/Al  junctions  have  been  performed  and  the  results 
compared  with  predictions  of  the  theory.  In  all  cases 
a  satisfactory  agreement  between  the  theory  and  the 
experiment  was  found.  The  inclusion  of  non-equilibrium 
phonons  in  metal  electrodes  has  been  suggested  as  a 
further  refinement  of  the  theory. 

Measurements  on  some  Mg/Mg  and  Mg/Au  junctions 
have  shown  an  additional  oscillatory  type  structure 
in  the  derivative  of  the  conductance.  The  occurrence 
of  the  structure  has  been  found  to  be  independent  of 
metal  electrodes,  and  barrier  layers.  Under  the 
influence  of  external  perturbations  such  as  warming 
the  junctions  to  room  temperature,  applying  DC  bias 
up  to  1  eV  and  passing  currents  along  the  metal  electrodes, 
the  properties  of  the  junctions  were  irreversibly  al¬ 
tered.  This  suggests  that  impurity  contamination  of 
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the  barrier  is  responsible  for  the  observed  structure. 

Tunneling  experiments  have  been  performed  in 
an  attempt  to  observe  the  quantization  of  allowed 
energy  levels  in  Bi  films.  The  structure  in  Bi/Bi 
j unctions , however ,  could  not  be  correlated  with  the 
film  thickness,  thus  rendering  the  results  negative. 

As  a  byproduct  of  the  investigation,  a  phonon  spectrum 
of  Bi  has  been  determined. 


vi 


ACKNOWLEDGEMENT  S 


I  wish  to  express  my  deep  gratitude  to  Dr.  J.G. 
Adler,  my  research  supervisor,  for  his  cheerful  encour¬ 
agement  and  guidance  throughout  the  course  of  this 
project.  My  introduction  into  the  world  of  an 
experimentalist  could  not  have  been  accomplished  without 
his  patience  and  instructive  help. 

I  wish  to  thank  Dr.  H.J.  Kreuzer  for  many 
illuminating  and  valuable  discussions  pertaining  to 
the  theory  of  zero  bias  anomaly.  I  am  also  thankful 
to  Dr.  P.S.  Turner  for  examination  of  films  by  an 
electron  microscope.  I  am  grateful  to  the  technical 
staff  of  the  Physics  Department  for  supplying  liquid 
helium  and  assisting  in  construction  of  evaporators 
and  cryostats.  My  special  thanks  are  due  to  Mr.  D.P. 
Mullin  for  his  help  and  cooperation  throughout  the 
course  of  my  work  in  the  laboratory.  I  also  wish 
to  thank  Mr.  R.  Teshima  for  his  computer  programming 
assistance . 

It  is  a  pleasure  to  acknowledge  the  alternating 
financial  support  of  the  Physics  Department,  National 
Research  Council  and  the  support  arising  from  the  Izaak 
Walton  Killam  Memorial  Fund. 

vi  i 


Finally,  I  am  thankful  to  Mrs.  Mary  Yiu  for 
her  patience  and  excellent  typing  performance  of  this 
thesis . 


viii 


TABLE  OF  CONTENTS 


Page 

CHAPTER  1  THEORETICAL  DISCUSSION  1 

1.1  Introduction  1 

1.2  Theory  of  elastic  tunneling  1 

1.3  Theory  of  inelastic  tunneling  10 

CHAPTER  2  THE  ZERO  BIAS  ANOMALY  14 

2.1  Introduction  14 

2.2  Non-equilibrium  theory  of  tunneling  18 

2.3  Description  of  the  model  19 

2.4  Scattering  mechanism  24 

2.5  Predictions  of  the  theory  22 

CHAPTER  3  EXPERIMENTAL  TECHNIQUES  50 

3.1  Sample  preparation  50 

3.2  Formation  of  tunnel  barriers  53 

3.3  Film  thickness  determination  56 

3.4  Temperature  measurement  and  control  58 

3.5  Production  of  low  temperatures  62 

3.6  Measurement  of  tunneling  charac-  62 

teristics 

CHAPTER  4  EXPERIMENTAL  RESULTS  69 

4.1  Introduction  69 

4.2  The  problem  of  the  background  69 

4.3  General  considerations  in  extracting  82 

tb,  t.  and  Tdrift. 

ix 


, 


Page 


CHAPTER  4 

(cont ' d) 

4.4 

Comparison  of  the  theory  with  the 
experiment 

86 

4.5 

Discussion 

120 

CHAPTER  5 

RESONANT  ELECTRON  TUNNELING 

128 

5.1 

Introduction 

128 

5.2 

Experimental  investigations  of 
oscillations 

131 

5.3 

Analysis  of  the  experimental  data 

159 

5.4 

Resonant  tunneling  through  impurity 
states 

164 

5.5 

Discussion 

172 

CHAPTER  6 

TUNNELING  INTO  BISMUTH 

179 

6.1 

Quantum  size  effect  in  Bi 

179 

6.2 

General  considerations 

184 

6.3 

Discussion  of  results 

185 

6.4 

Phonons  in  Bi 

193 

6.5 

Other  interesting  effects 

202 

CHAPTER  7 

SUMMARY  AND  CONCLUSIONS 

206 

7.1 

Discussions 

206 

7.2 

Suggestions  for  further  work 

209 

BIBLIOGRAPHY 

212 

APPENDIX  I 

219 

APPENDIX  II 

226 

x 


. 


LIST  OF  TABLES 


Page 

Table 

i — 1 

• 

CO 

Selection  of  materials  for 
evaporation 

54 

Table 

1 — 1 

• 

Thickness  and  periodicity  of 
oscillations 

182 

Table 

6.2 

Identification  of  phonon  peaks 
in  a  Bi/Bi  junction 

200 

Table 

6.3 

Identification  of  phonon  peaks 
in  Mg/Bi 

201 

Table 

6.4 

Identification  of  phonon  peaks 
in  Al/Bi 

201 

x  l 


LIST  OF  FIGURES 


Figure 

Page 

1.1 

Transfer  Hamiltonian  model. 

3 

1.2 

Schematic  energy  level  diagram  for 
tunneling  between  normal  metals. 

11 

2.1 

Conductance  curves  for  junctions  with 
various  electrodes. 

15 

2.2 

Schematic  diagram  of  scattering 
processes . 

20 

2.3 

Illustration  of  a  deblocking  by  mag¬ 
netic  field. 

28 

2.4 

6a  vs  V  and  —  vs  V  at  T  =  0  K 

(a2F  =  alU2).°°  dV 

36 

2.5 

S,a2lS- V  aD  dV  vs  V  at  T  0  K 

(cx  F  a  2  GJ  ). 

38 

2.6 

2 

6a  vs  T  at  V  =  0  mV  (a  F  =  a^aO  . 

42 

2.7 

2.8 

2.9 

2  2 

6a  vs  T  at  V  =  0  mV  (a  F  =  a ) . 

i 

Vm-,„  vs  T  and  - - -srr-  vs  T 

max  On  dV  v 

/  2^  ,  °  vmax 

(a  F  =  a-,  oo)  . 

1  d6oe 

V  vs  T  and  —  — ttt—  vs  T 

max  o0  dV  vmax 

(a  F  =  a2oo  )  . 

44 

46 

48 

3.1 

Layout  of  the  sample. 

51 

3.2 

Schematic  diagram  of  the  cryostat. 

60 

3.3 

Simplified  bridge  configuration. 

63 

3.4 

Schematic  diagram  of  the  data 
acquisition  system. 

67 

4.1 

a  vs  V  for  a  Mg/Mg  junction  for  various 
temperatures . 

70 

4.2 

~~  ~  vs  V  for  a  Mg/Mg  junction 
aQ  dV 

for  various  temperatures. 

72 

xi  1 


. 


' 


Figure 


Page 


4.3 

4.4 

4.5 

4.6 

4.7 

4.8 

4.9 

4.10 

4.11 


4.12 

4.13 

4.14 


6a  vs  V  for  a  Mg/Mg  junction  for 
various  temperatures. 


1  d6a 
ao  dV 
various 


vs  V  for  a  Mg/Mg  junction  for 
temperatures . 


Experimental  a0^/G0  vs  V' 


Theoretical  fit  to  experimentally 

1  d6ae 

determined  Vmax  vs  T  and  — 


vs  T  for  a  Mg/Mg  junction. 


o 


dV 


V. 


max 


75 

77 

80 

87 


Theoretical  fit  to  experimentally 
determined  6a/aQ  vs  T  for  a  Mg/Mg 
junction  at  V  =  0  mV. 


Theoretical  and  experimental 

d6a 

■ —  vs  V  for  a  Mg/Mg  junction. 

g0  dV 

a  vs  V  for  085A  Mg/N2/Au  for  various 
temperatures . 

a  vs  V  for  085B  Mg/N^/Au  for  various 
temperatures. 

Theoretical  fit  to  experimentally 

determined  Vraax  vs  T  and 

d5a.. . 

vs  T  for  085B1  and 


a 


dV 


o 

085B2. 


V. 


max 


Theoretical  and  experimental 
,  d6a 

— - vs  V  for  a  Mg/N0/Au  junction 

oQ  dV  2 

So/oQ  vs  T  at  V  =  0  mV,  for  two 
junctions  with  different  thickness 
of  metal  electrodes. 

Theoretical  fit  to  experimentally 
determined  Vm_v  vs  T  and 

aCa-  vs  T  for  an  Al/Al 


1 


dV 


V 


max 


junction . 


90 


92 


95 

97 

99 


101 


104 


107 


xiii 


Figure 


Page 


4.15 

6g  vs  V  for  an  Al/Al  junction  for 
various  temperatures. 

109 

4.16 

Theoretical  fit  to  experimentally 
determined  Vmax  vs  T  and 

1  d5Gp 

—  - ■  T7  vs  T  for  an  Al/Al 

°o  dv  vmax 

junction . 

112 

4.17 

Theoretical  and  experimental 
i  d6aP 

—  ■  •  vs  V  for  an  Al/Al 

0o  dV 

junction . 

114 

4.18 

Illustration  of  the  appearance  of 
the  energy  gap  in  an  Al/Al  junction. 

116 

4.19 

a  (0 )  vs  T  for  851  Al/Al,  691  Al/Al. 

118 

5.1 

—  vs  V  for  877A  Mg/Au. 

a  dV 
o 

129 

5.2 

a  vs  V  for  883A  Mg/Mg. 

132 

5.3 

i|!vs  V  for  883A  Mg/Mg . 

0o  dV 

134 

5.4 

Temperature  smearing  of  oscillations 
in  867B  Mg/Au. 

138 

5.5 

o  vs  V  for  875B  Mg/Au. 

140 

5.6 

o/o  vs  V,  lda/aQdV  vs  V  for  875  Mg/ 
Mg/Su . 

142 

5.7 

a  vs  V  for  95A  Mg/^/Au. 

145 

5.8 

_L  4?  vs  V  for  95A  Mg/N_/Au. 
aQ  dV  ^ 

147 

5.9 

a/a  vs  V  for  080B  Mg/Au  before  and 
after  annealing. 

151 

5.10 

-i-  vs  V  for  080B  Mg/Au  before 

aQ  dV 

and  after  annealing. 

153 

5.11 

a/a  (0)  vs  V  for  J107B  Mg/N2/Au  before 
and  after  applying  high  bias  across 
the  junction. 

155 

xiv 


' 


Figure 


Page 


5.12 

—  ~  vs  V  for  J107B  Mg/No/Au 
aG  dV  ^ 

before  and  after  applying  high 

bias  across  the  junction. 

157 

5.13 

Schematic  potential  diagram  for 
a  junction  with  an  impurity  lying 
inside  the  barrier. 

165 

5.14 

oR/oN  vs  V. 

169 

5.15 

Spin  splitting  of  an  impurity 
level . 

174 

6.1 

o/oQ  vs  V  for  001,  007,  016  Bi/Bi 
junctions . 

186 

6.2 

i  ^  vs  V  for  001,  007,  016  Bi/Bi 

0o  dV 

junctions . 

188 

6.3 

~  vs  V  for  007  Bi/Bi. 

oG  dv 

191 

6.4 

Phonons  in  Bi/Bi. 

194 

6.5 

Phonons  in  Mg/Bi. 

196 

6.6 

Phonons  in  Al/Bi. 

198 

6.7 

Normalized  conductance  of  typical 
Al/Bi,  Mg/Bi  and  Bi/Bi  junctions. 

203 

xv 


CHAPTER  1 


THEORETICAL  DISCUSSION 

1 . 1  Introduction 

The  concept  of  a  particle  tunneling  from  a  classi¬ 
cally  allowed  region  of  space  through  a  classically 
forbidden  region,  into  the  second  allowed  region  was  one 
of  the  early  predictions  of  Schrodinger ' s  wave  equation. 

The  possibility  that  a  tunnel  current  may  flow 
between  metals  separated  by  a  vacuum  was  first  considered 
by  Frenkel  (Frenkel,  1930).  Sommerfeld  and  Bethe 
(Sommerfeld  and  Bethe,  1933)  extended  Frenkel's  calcu¬ 
lation  to  the  case  in  which  the  vacuum  was  replaced  by 
a  thin  insulating  layer.  The  actual  physical  fabrication 
of  such  a  system  was  delayed  until  1960  when  the  thin 
film  technology  had  sufficiently  advanced  allowing  Fisher 
and  Giaver  (Fisher,  Giaver,  1961)  to  build  a  tunnel 
structure  using  aluminum-aluminum-oxide-aluminum.  The 
results  of  their  experiments  stimulated  a  sequence  of 
measurements  on  both  normal  and  superconducting  tunnel 
junctions  with  far  reaching  results  leading  to  a  better 
understanding  of  electron  transport  phenomena  in  metals. 

1 . 2  Theory  of  elastic  tunneling 

Bardeen  (Bardeen,  1961)  has  regarded  the  tunneling 


1 


, 

■ 
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process  as  a  transition  of  an  electron  across  the 
barrier  which  itself  acts  as  a  perturbing  potential. 

The  composite  metal-insulator-metal  system  is  described 
by  Hamiltonians  H^,  of  the  left  and  right  hand 
metals  and  one  electron  eigenstates  ^  and  ip2  which 
are  weakly  coupled  together  by  a  perturbation  U  (fig. 1.1). 
By  taking  a  linear  combination  of  ip^  and  ip2  we  seek  then 
a  modified  solution  to  the  time  dependent  Schrodinger 
equation.  Knowing  i p^  and  ip2  the  transition  matrix 

element  °f  the  Hamiltonian  between  states  ip^  and  2 

can  then  be  calculated  using  the  method  of  time  depen¬ 
dent  perturbation  theory. 

The  transition  probability  P^2  per  unit  time  of 
an  electron  in  a  state  to  a  state  \p2  is  then  given  by 
Fermi  golden  rule 


P 


12 


2tt 

fi 


p2f°d-f°) 


d.i) 


where  p 2  is  the  density  of  final  states  in  the  metal  2, 
f°  and  f°  are  the  Fermi-Dirac  probabilities  of  occupa- 

_L  *—> 

tion  of  states  ip^  and  ip2> 

As  an  approximation  for  the  electron  current 
density  between  the  two  metals,  an  independent  particle 
model  is  used.  It  is  assumed  that  the  metal-insulator 
interface  is  planar  and  that  the  component  of  the  wave 
number  parallel  to  the  plane  of  the  barrier ,  k,  , ,  is 
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Fig.  1.1 


Transfer  Hamiltonian  Model 
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conserved  before  and  after  the  transmission.  The  total 
energy  of  the  electron  of  course, is  conserved  and  this 
is  implied  when  writing  eq.  (1.1).  To  obtain  the  tunnel¬ 
ing  current  from  metal  1  to  metal  2'  P12  1S  summed 
over  all  states  1  of  fixed  ,  summed  over  k^  ,  multi¬ 
plied  by  2  for  spin  and  by  charge  e  of  the  electron. 

The  sum  over  states  1  may  be  replaced  by  an  integral 
over  energy  at  fixed  wave  number  .  Subtracting  the 
current  from  metal  2  to  metal  1 ,  one  obtains  for  the  net 

oo 

’  |M12|2Pl(Ex)p2(Ex+eV)  [f°(E)-f°(E+eV)]dEx 
0 

(1.2) 


current 
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where 


f 


o 


(E) 


1 

(E-Ef  )/K  T 

1  +  e 


and  p.,  i  =  1,2,  represents  a  one  dimensional  density 
of  states,  associated  with  the  motion  perpendicular  to 
the  barrier 


pi  = 


L  ,  9E  > 

^  3kx. 

1 


-1 


1,2 


(1.3) 


where  L  gives  the  length  of  the  metal.  Unless  otherwise 
stated  all  energies  are  measured  with  reference  to  the 
bottom  of  the  conduction  band  of  metal  1.  Because  of 
the  independent  particle  model,  we  can  write  E  as 
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E  =  E  +  E 
x 


Wlth  *2,2  *2.2 
■h  k  ft  k 

p  —  _ Li  p  =  - - 

h\\  2m  '  x  2m 

i  1 2 

The  functional  form  of  the  matrix  element  |M^2  I 
depends  on  the  profile  of  the  barrier  potential  used 
for  the  calculation  of  electron  transmission  probabi¬ 
lity.  Harison ,  (Harrison,  1961)  applying  a  WKB 
approximation  found  that 


D(E,Ej  ) 
7“2 

4tt  Plp2 


(1.4) 


where  the  barrier  transmission  coefficient  D  is  given 
by 


D (E,E,  , )  =  exp  (-2 


0 


~  (U(x)-(E-E  )dx) 

n  1 1 


(1.5) 


2 

Substituting  eq .  (1.4)  for  |M12|  into  eq.  (1.2)  we 

obtain 
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2e 

-h 


i 
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D(E,En  )  (f°(E)  -  f°(E+  eV)  dE 

0 


(1.6) 


and  find  that  the  density  of  states  factors  appearing 
in  eq.  (1.2)  have  cancelled  with  the  matrix  element  so 
that  the  tunneling  current  does  not  explicitly  depend 
on  them.  The  absence  of  density  of  states  in  eq .  (1.6) 

is  a  result  of  the  independent  particle  model,  which 
enabled  the  separation  of  the  electron  wave  number  into 


parallel  and  perpendicular  components  with  respect  to 
the  plane  of  the  barrier  and  the  replacement  of  the 
sum  over  the  normal  components  by  an  integral  over 
energy  E  .  Furthermore,  the  WKB  approximation  is 

restricted  only  to  regions  where  the  fractional  change 

(5  A 

in  the  electron  wavelength  is  small  i.e.  -j-  <<  1. 

This  is  not  satisfied  for  electrons  near  the  band  edge, 
having  small  kinetic  energy.  Therefore,  as  Harrison 
points  out,  band  edges  and  energy  gaps  may  introduce 
an  observable  structure  in  the  tunneling  current. 

The  form  of  the  matrix  element  eq.  (1.3)  and 
its  use  in  eq .  (1.5)  implies  that  the  tunneling  is 

from  one  quasi-continuous  band  of  states,  through  the 
barrier  into  another  quasi-continuous  band.  If,  on 
the  other  hand,  p2  is  a  set  of  localized  surface  states 
or  impurity  states,  the  resulting  matrix  element  is  no 
longer  proportional  to  p2^  but  rather  to  t  ^ ,  the  life¬ 
time  of  the  state  (Penley,  1962) ,  and  the  density  of 
such  states  does  not  cancel.  In  general,  if  the  rate 
limiting  quantity  in  the  tunneling  process  is  the 
electron  velocity,  and  if  the  velocity  is  inversely 
proportional  to  the  density  of  states,  which  is 
characteristic  of  a  band  of  single  particle  states, 
then  the  tunneling  current  does  not  contain  information 
about  the  density  of  states.  But,  if  the  tunneling 
process  is  primarily  via  traps  and  defect  states  in 
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the  insulator,  which  can  be  characterized  by  a  life¬ 
time  t,  the  density  of  states  should  not  cancel. 

Eq.  (1.6)  can  be  cast  into  a  more  convenient 

form  by  converting  the  sum  over  k . ,  to  an  integral 

I  I 

over  k||  and  then  to  an  integral  over  E||  ,  with 

d2k  M  =  dE,  . 

I!  27Tfl  " 
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oo 


47TmeA 
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[f°(E) 


f  °  (E  +  eV)  ]  D  (E  ,E  . ,  )  dE  dE 

I  I  X 


(1.7) 


The  upper  limit  on  the  integration  over  dE | ^  can  be 
determined  from  the  condition  that  the  component  of 
the  wave  number  parallel  to  the  plane  of  the  barrier 
k  |  and  the  electron  energy  E  are  conserved  during 
transition.  Exchanging  the  order  of  integration  in 


eq . 

(1.7) 

and  noting  that  D  (E , E |  | 

)  =  D  (E  )  and 

.A. 

3E  ! 

,  /9E  = 

1  we  obtain 

00  E 

I  = 

4TTmeA 

>3 

[f°  (E)  -  f°  (E  +  eV)  ]  dE 

D  (E  )  dE 

X  X 

(1.8) 

where  A  is  the  effective  junction  area. 

For  illustration,  we  calculate  the  tunneling 
current  at  T  =  0  K.  Since 


f° (E)  -  f°  (E  +  eV)  =0  , 


unless 


eV  <  E  <  Efl 


when 


' 
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f° (E)  -  fu  (E  +  eV)  =  1 


■  o 


eq.  (1.8)  reduces  to: 


I  = 


4TrmeA 

h3 


Ef i~eV 


dE 


E 


0 


D (E  )  dE 
x  x 


(1.9) 


or 


E 


I=  IrrmeA  (eV) 


rfl 


D (E  )  dE 
X 


(1.10) 


0 


providing  that  eV  <<  E^^.  The  current  is  therefore 
ohmic  at  low  biases.  The  current  at  T  ^  0  K  can  be 
connected  to  the  current  at  T  =  0  K  (Stratton,  1962) 


by 


I  (V,T) 
I (V,0) 


*ckbt 


sin  ttCKbT 


(1.11) 


where  C  is  a  constant  dependent  on  the  barrier  height 
and  thickness.  The  expected  variation  of  the  current 
with  temperature  is  more  apparent  in  eq.  (1.11) 
(Maissel ,  Glang,  1970) 


-9  ( £ . T ) 


$ 


(1.12) 


where  £  is  in  A,  and  T  and  $  (average  barrier  potential) 
are  in  conventional  units.  For  typical  parameters 
(j)  =  1  eV  ,  £  =  30  A,  T  =  5  K, 

AI(V,T)  %  10"5 


which  is  negligible. 


. 
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1 . 3  Inelastic  tunneling 

Tunneling  which  involved  a  loss  of  energy  by 
tunneling  to  localized  impurities  in  the  barrier  was 
first  reported  by  Jacklevic  and  Lambe  (Jacklevic  and 
Lambe,  1966).  In  studying  Al-insulator-Pb  tunnel 
junctions,  they  have  noticed  increases  in  the 
conductance  of  the  junctions  at  characteristic  voltages 
corresponding  to  vibrational  frequencies  fia)  of  molecular 
species  contained  in  the  insulating  region  of  the 
junction.  The  energy  level  diagram  on  fig.  1.2  allows 
a  visualization  as  to  how  such  processes  can  lead  to 
increases  in  conductance. 

Fig.  1.2  represents  a  tunnel  junction  energy 
level  diagram  where  the  Fermi  levels  of  the  two  metals 
are  separated  by  the  energy  eV  where  V  is  the  applied 
voltage.  For  simplicity  let  us  assume  that  the  system 
is  at  T  =  0  K  which  means  that  the  Fermi  levels  are 
sharp  and  the  impurity  present  in  the  barrier  is  in 
its  ground  state.  Electrons  may  tunnel  elastically 
from  metal  1  to  metal  2  (process  a) .  If  a  new  mechanism 
exists  so  that  an  electron  can  lose  energy  ftm  to  a 
localized  impurity  in  the  barrier  as  it  is  tunneling 
(process  b)  a  new  channel  will  open  up.  This  new 
channel  is  possible  only  if  there  exist  empty  states 
in  metal  2  at  an  energy  -hoo  less  than  This  becomes 

possible  when  eV  >  ii m . 


' 


Fig.  1.2 


Schematic  energy  level  diagram  for 
tunneling  between  normal  metals: 

a)  elastic  tunneling 

b)  inelastic  tunneling. 
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This  onset  of  a  new  tunneling  channel  results  in 
a  steplike  increase  in  the  conductance  and  6  function 
like  structure  in  the  derivative  of  the  conductance. 

A  finite  temperature  will  smear  the  Fermi  surface  of 
both  metals  and  thus  it  will  smear  the  6  function  like 
structure  although  it  will  be  centered  around  eV  =  hw . 
The  energy  loss  of  the  tunneling  electron  is  not  limited 
only  to  losses  to  localized  impurities.  In  fact,  as 
it  was  shown  by  Giaver  and  Zeller  (Giaver  and  Zeller, 

1968)  and  others,  inelastic  excitations  due  to  phonons 
of  the  barrier  (Adler,  1969),  phonons  of  metal  electrode 

(Rowell  et  al,  1969)  itself,  surface  plasmons  (Tsui, 

1969) ,  magnons  (Tsui,  et  al,  1971)  are  possible. 


CHAPTER  2 


THE  ZERO  BIAS  ANOMALY 

2 . 1  Introduction 

The  description  of  a  tunnel  junction  at  small 
biases  (eV  <<  $)  as  an  ohmic  element  (eq.  1.10)  is 
found  to  be  adequate  only  when  analyzing  I-V  curves. 
Conductance  measurements  using  derivative  techniques 
to  an  accuracy  of  ^1  %  have  shown  that  the  conductance 
of  a  tunnel  junction  up  to  500  meV  is  roughly  para¬ 
bolic 

a (V)  =  a  +  +  yV^  .  (2.1) 

Theoretical  calculations  (Brinkman  et  al ,  1969)  have 
explained  the  offset  of  the  parabola  from  zero  bias 
as  being  due  to  the  barrier  asymmetry  arising  from 
different  work  functions  of  the  metal  electrodes  . 
Offsets  larger  than  50  meV'  were  attributed  to  the 
distortion  of  the  barrier  potential  due  to  impurities 
trapped  in  the  junction.  The  parabolic  fit  by  eq. 
(2.1)  is  considered  to  be  approximate  in  the  sense 
that  it  helps  to  determine  the  background  conductance 
level  for  any  additional  structure  which  may  appear 
superimposed  on  the  tunnel  characteristics.  Fig.  2.1 
illustrates  such  a  case  (Chen,  Adler,  1970).  The  set  of 
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Fig .  2.1 


Normalized  condutance  curves  for  junctions  with 
various  electrodes.  The  curves  are  displaced 
vertically  by  0.1  for  clarity.  The  scales  for 
Al-Al,  Al-Mg  junctions  should  be  multiplied  by 
two.  The  insert  shows  the  sizes  of  zero  bias 


anomalies  of  some  junctions. 
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curves  for  different  junctions  is  roughly  parabolic 
with  weak  conductance  increases  at  118  meV  due  to 
inelastic  excitation  of  aluminum  hydrate  (Geiger  et 
al ,  1969). 

The  insert  in  fig.  2.1  shows  a  structure  around 
zero  bias  (commonly  called  the  zero  bias  anomaly) ,  the 
detection  of  which  reguires  conductance  measurements 
with  a  resolution  to  one  part  in  104.  Disregarding 
the  relatively  large  conductance  dip  for  the  Al/Ni 
junction ,  which  is  associated  with  the  magnetic  pro¬ 
perties  of  the  nickel  electrode  (Rowell,  1969),  small 
conductance  decreases  (<  1  %)  around  zero  bias  in  normal 
metal-insulator- junctions  (seen  here  in  Al/Al  and  Al/Ag 
junctions)  were  for  a  long  time  unexplained.  The 
structure  does  not  fall  into  any  other  previously 
investigated  classes  of  zero  bias  anomalies  explained 
on  the  basis  of  metallic  inclusions  (Giaver ,  Zeller, 
1968)  or  magnetic  impurities  inside  the  insulating 

layer  (Rowell,  Shen,  1967). 

A  recent  theory  (Trof imenkof f  et  al ,  1972) 
through  a  nonequilibrium  treatment  of  the  tunneling 
process  has  shown  that  the  appearance  of  the  low  bias 
conductance  minimum  in  M-I-M  junctions  is  closely 
related  to  the  finite  electron  relaxation  rates  in 


metal  electrodes. 
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In  what  follows,  we  shall  present  the  theoretical 
model  of  describing  the  non-equilibrium  tunneling  process 
with  its  predictions  which  are  relevant  to  experimental 

verification . 

2 . 2  Non-equilibrium  theory  of  tunneling 

A  conventional  theory  of  tunneling  (Chapter  1) 
assumes  that  the  tunneling  electron  is  in  thermal 
equilibrium  with  the  rest  of  the  electrons  inside  the 
metal.  This  is  expressed  by  using  the  equilibrium 
Fermi-Dirac  distribution  function,  f ° ,  in  the  calcula¬ 
tion  of  the  tunneling  current. 

Consider,  however,  the  following  argument.  An 

electron  which  penetrates  the  barrier  finds  itself  at 

4 

certain  energy  level  in  the  host  metal.  At  He  temper 
atures,  the  final  state  of  the  electron  is  such  that 
even  at  small  imposed  biases  ^1  meV ,  it  is  energetically 
well  above  the  thermally  smeared  Fermi  sea  (1  K = 86  pV) . 
The  existing  non— equilibrium  situation  is  neutralized 
only  after  a  short  period  of  time  which  is  character 
ized  by  the  relaxation  processes  in  the  metal.  During 
the  time  the  excited  state  is  occupied,  by  the  virtue 
of  Pauli  exclusion  principle,  all  incoming  electrons 
are  effectively  blocked  from  entering  the  same  state. 

If  the  decay  times  of  the  excitation  were  of  the  same 
order  of  magnitude  as  the  time  an  electron  decides 


, 
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to  penetrate  the  barrier,  the  blocking  of  the  available 
states  would  manifest  itself  by  a  strong  decrease  in  the 
conductance  of  the  tunnel  junction.  The  use  of 
equilibrium  Fermi-Dirac  distribution  function  is  then 
equivalent  to  an  assumption  that  the  thermalization  is 
immediate.  In  other  words,  the  relaxation  times  of  the 
interactions  are  assumed  to  be  equal  to  zero. 

To  obtain  quantitative  predictions  of  the  above 
mentioned  process,  a  non- equilibrium  distribution  func¬ 
tion  has  to  be  found,  characteristic  of  the  existing 
situation  in  the  metal  electrode.  This  is  done  by 
solving  a  transport  model  for  the  tunnel  current,  based 
upon  the  Boltzmann  equation. 

2 . 3  Description  of  the  model 

With  reference  to  fig.  2.2,  assume  that  an 
electron  tunnels  from  a  state  ij^  into  a  state  \p2-  By 
the  virtue  of  Pauli  exclusion  principle  the  final 
state  1^2  is  unavailable  for  tunneling  as  long  as  it 
is  occupied.  However,  through  various  scattering 
processes,  the  electron  is  scattered  out  of  the  occu¬ 
pied  state.  The  competing  mechanisms  of  filling  and 
emptying  the  states  soon  establish  across  the  junc¬ 
tion  a  steady  state  in  which  the  rate  of  change  of 
the  distribution  function  with  time  is  zero,3f/Bt  =  0. 


Fig .  2.2 

Schematic  diagram  of  scattering  processes: 

a)  electron  in  a  forward  direction 

b)  inelastic  scattering  by  phonon  emission 
or  absorption 

c)  ,  d)  elastic  impurity  scattering 
e)  specular  boundary  scattering 
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Balancing  the  rate  of  change  of  the  distribution  func¬ 
tion  in  metal  2  due  to  incoming  (8f2/3t)1^2  and  outgoing 
( 3f 2/8t) 2^^)  electron  flux  with  losses  due  to  scatter¬ 
ing  processes,  (3f2/3t ) SCatt  we  ^ave: 


9  f . 


3 1 


-)  =  0  =  ( 


3f0  3fo 

+  (3ti)  +  (3t-> 

3t  i->2  2+1  scatt 


(2.2) 


In  relaxation  time  approximation 


(ii) 

'  3 1 ' 


scatt 


Af 

T 


(2.3) 


where.  Af  is  the  deviation  of  f  from  equilibrium 


f  =  f°  +  Af 


(2.4) 


and  x  is  a  time  characteristic  of  scattering  processes 
which  bring  f  back  to  equilibrium.  With  the  help  of 
Fermi  golden  rule,eq.  (2.2)  is  written  as 


3  f . 
3 1 


2  TT 

h 


M 


12 


Plfl(l  f2^ 


2tt 

h 


M 


12 


plf 2  ^1 


fl} 


Af. 


(2.5) 


or 


Af- 


2  TT 


M. 


T 


1 T  I  “12  I  plv"l  J"  2 


(fn  -  fo) 


(2.6) 


By  symmetry  a  similar  equation  is  written  for  f^ 


Af 


1 


T 


1 


2  IT 

h 


p2  2 


(2.7) 
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Since  the  relevant  quantity  which  enters  the  calcula¬ 
tion  of  the  tunneling  current  is  the  difference  in  the 
distribution  functions  (eq.  1.8),  the  two  equations 
above  are  solved  in  that  form: 


fl  f  2 


-O 

1 


-O 


(2.8) 


n  2  7T  | 

1  +  TdM12 


[  Pf1 2  P  2  T1 


On  the  basis  of  experimental  evidence,  we  expect  that 
the  deviations  from  equilibrium  are  small.  This  permits 
to  rewrite  eq.  (2.8) 


fl  f  2 


(f°  -  f°)  (1 


2  TT 


M. 


h  '“12 


[  P 1^  2 +  P2T1^  ^  ' 


(2.9) 


Therefore  instead  of  the  usual  difference  in  the 
equilibrium  distribution  functions  (f°-  f^)  •  (ec3*  d*8), 

the  above  relationship  between  f^  and  should  be  used 
in  the  calculation  of  the  tunneling  current.  In  the 
next  section  we  shall  consider  the  nature  of  the 
scattering  processes,  in  order  to  determine  the  different 
relaxation  mechanisms  present  in  the  metal. 


. 


2  4 


2 . 4  Scattering  mechanism 

An  electron  in  the  metal  undergoes  a  number  of 
collisions  which  in  general  may  alter  its  energy  and 
momentum.  Inelastic  scattering  due  to  phonon  emission 
and  absorption  is  characterized  by  a  temperature  and 
energy  dependent  relaxation  time  t(£,T) .  A  first  order 
calculation  applied  to  an  interacting  electron-phonon  gas 
(Wilkins ,  1968)  gives  the  following  result: 


1 


(e,T) 


d  (hoo)  a* 2  (ha))  F  (hco)  [2N°  (had  + 

4 

0 


+  l+f°(e+hto)  -  f0(e-fiw)] 


(2.10) 


where  e  is  the  excitation  energy  measured  relative  to 

the  Fermi  level.  a  F  is  a  product  function  involving 

2 

the  electron-phonon  coupling  function  a  (ftw)  and  the 
phonon  density  of  states  F  (haj)  .  The  equilibrium  dis¬ 
tribution  function  for  the  phonons  is  given  by 


ekfiT  - 

Since  the  evaporated  films  are  polycrystalline  and 
thus  have  a  characteristic  of  an  isotropic  material, 
the  relaxation  time  is  assumed  to  be  isotropic. 

To  facilitate  a  comparison  of  the  theory  with 
the  experiment  it  is  necessary  to  investigate  the 


25 


behaviour  of  t  at  low  temperatures  and  low  excita- 

ep 

2 

tion  energies.  Approximating  a  F  at  low  frequencies 
by  a  power  law  behaviour 


a 


anm 


n 


n  =  integer 

a  =  constant 
n 


Wattamaniuk  (Wattamaniuk ,  1972)  found  that 


i  ^a  -j  j  2 

1 ,  -r  =  (e  +  (ttK  T)  ] 

Tep  £  B 


for  n  =  1,  whereas  for  n  =  2 


,  2-jia9 

— 7^(3)  (kRT) 

( e )  _^n  ft  B 

ep 


1 


Tep(e)  T+0 


27Ta2  H3 

~W~  3 


(2.12) 


(2.13) 


(2.14) 


where  E,  (x)  is  a  Riemann  zeta  function. 

The  two  limiting  cases  of  x  appear  as  a  result 

of  a  controversy  which  at  the  present  time  exists  on 

the  low  energy  limit  of  a  F.  At  low  frequencies  F(d) 

2 

in  the  Debye  approximation  varies  as  w  •  Therefore 
for  n  =  2,a^(fLoo)  is  a  constant.  Experimental  results 
(Goy,  Castaing, 1973)  point  to  this  behaviour. 
Theoretical  calculations  (Carbotte,  Dynes,  1968;  Allen, 
Cohen,  1970)  indicate  a  possible  linear  behaviour  of 
a^F  due  to  Umklapp  processes.  We  shall  consider  both 


. 


cases  separately,  hoping  that  the  experimental  results 
will  favour  one  limit. 
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Lacking  detailed  knowledge  about  the  structure 
of  evaporated  films,  scatterings  due  to  lattice  defects, 
impurities  and  diffuse  boundaries  are  characterized  by 
a  constant  relaxation  time  x^  i.e.  the  scattering  is 
considered  to  be  energy  and  temperature  dependent. 

Another  mechanism  by  which  the  blocking  of  the 
incoming  electron  flow  is  decreased, is  obtained  by  an 
application  of  a  magnetic  field  perpendicular  to  the 
direction  of  the  tunnel  current. 

An  electron  which  penetrates  the  barrier  will 
travel  away  from  it  and  can  propagate  across  the  film 
before  being  scattered.  Assuming  specular  reflection 
at  the  boundary,  the  electron  is  reflected  back  with  a 
reversed  momentum.  Upon  reaching  the  metal-insulator 
interface,  it  may  again  be  reflected  specularly,  thus 
reversing  the  momentum  back  to  the  original  one.  If 
the  travel  time  over  the  length  of  the  film  is  much 
less  than  the  time  an  electron  decides  to  penetrate 
the  barrier  or  the  scattering  time,  then  the  electron 
effectively  remains  in  its  quantum  state  and  contributes 
to  the  blocking.  For  typical  metals  at  energies  close 
to  the  Fermi  level  Vf ^  108  cm/sec  and  for  a  1000  K 
film,  the  travelling  time  back  and  forth  is  ^2x10 


sec 


■ 


■ 


1 
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which  is  at  least  two  orders  of  magnitude  less  than 

the  relaxation  times  due  to  electron-phonon  inter- 
-11  4 

action  (^10  sec  at  He  temperatures). 

In  the  presence  of  the  magnetic  field  (Fig.  2.3)  , 
oriented  perpendicularly  to  the  direction  of  the 
tunnel  current  flow,  the  Lorentz  force  induced  drift 
provides  a  mechanism,  whereby  the  specularly  reflected 
electrons  move  down  the  electrode  and  ultimately  move 
away  from  the  tunneling  area.  In  a  classical  picture, 
the  relaxation  time  associated  with  this  process  is 
then  equated  to  the  time  it  takes  an  electron  to  drift 
out  of  the  junction  area  and  thus  be  removed  from  the 
forward  blocking  direction.  It  is  assumed  that  at  low 
biases  the  magnetic  field  introduces  no  appreciable 
changes  in  the  Fermi  surface  of  the  metal  which  could 
affect  the  tunnel  current.  The  validity  of  this 
assumption  is  confirmed  by  the  experimental  results. 

With  reference  to  Fig.  2.3,AA  represents  the 
metal-insulator  interface,  which  is  at  distance  L 
away  from  the  other  end  of  the  film.  A  magnetic  field 
of  B  kG  produces  a  circular  motion  with  radius 

mcVr 

_  r  t  o  i  c.  \ 


where  c  is  speed  of  a  light,  and  is  the  velocity  of 
the  electron  at  the  Fermi  level.  During  the  time 


' 


Fig .  2.3 

Schematic  illustration  of  magnetic 
field  induced  drift  velocity  of  a 
tunneling  electron. 
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T  =  —  ,  the  electron  drifts  down  the  electrode  a 
L  AB  Vf 

distance  S 


S 


(2.16) 


Approximating  AB  by  a  chord  AB 

AA  =  /s2  + 

becomes: 

AB 

_  =  /s2  +  L2 

TAA  Vf 


(2.17) 


(2.18) 


The  drift  velocity  is  given  as: 

SV. 


V 


drift 


T 


AA 


S2  +  L2 


(2.19) 


and  thus  if  w  denotes  the  width  of  the  film; 

7  -1/2 


-  /Vdrif  t  - 
T drift  w 


w/2 
Vf  l 


1  -A  -  K 

R  • 


(2.20) 


For  L  =  1000  A  and  R(60  kG)  =  20000  A,  L/R  <<  1  and 
we  can  expand  the  square  root  in  eq.  (2.20)  in  powers 
of  L/R.  Keeping  the  first  non-zero  term  we  have: 


Tdrif t (B) 


2mc  ,w> 
eB  L 


(2.21) 


The  final  result  shows  that  the  effect  of  a  magnetic 
field  is  to  reduce  the  blocking  by  a  rate  which  depends 
entirely  on  the  geometry  of  the  junction  via  the  (j-j)  ratio. 


, 


A  magnetic  field  applied  parallel  to  the  tunnel 
current  provides  no  Lorentz  force  on  the  electron  and 
thus  should  have  no  effect  on  the  size  of  the  zero  bias 
anomaly . 

Recapitulating,  we  have  considered  three  possi¬ 
ble  mechanisms  which  may  scatter  electrons  from  a 
given  state: 

a)  Energy  and  temperature  dependent  electron-phonon 
scattering  time. 

b)  Temperature  and  energy  independent  impurity 
scattering  time,  which  includes  contributions 
from  lattice  defects  and  diffuse  boundaries. 

c)  Scattering  due  to  magnetic  field  induced  drift 
velocity  of  the  tunneling  electron. 

Since  the  probability  of  a  scattering  is  inverse 
proportional  to  the  relaxation  time  and  since  the  dif¬ 
ferent  scattering  mechanisms  are  independent  of  each 
other,  we  can  write 


1 
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E-  +  A_  + 


ep 


T  . 
1 


T 


drift 


(2.22) 


where  x  is  the  effective  relaxation  time  of  an  excited 


electron . 


. 


32 


2 . 5  Predictions  of  the  theory 

Substituting  eq.  (2.9)  into  eq .  (1.8)  the  total 

tunneling  current  I  is  given  as 


I  =  I  +  61  , 
o 


(2.23) 


where  I  is  the  equilibrium  current  and  61  is  the 
o 

blocking  current 
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4TrmeA 


(f° (E)  -  f° (E  +  eV) )  dE 


E 

x 


'  °2(V 

✓ 
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f  1  T1^E-Efl^  1 2  (E-Efl+eV)« 
2-rrh  p  1  (Ex )  +  2iTh  p2(Ex  +  eV) 


(2.24) 


When  writing  the  above  equation  we  assume  that  the 
energy  dependence  of  the  relaxation  time  is  coming  only 
from  the  energy  dependence  of  the  electron-phonon 
counterpart . 

Notice  the  dependence  of  6 1  in  the  square  of 
the  barrier  transmission  coefficient,  thus  making  61 
a  sensitive  function  of  the  insulator  thickness. 
Furthermore  the  density  of  states  factors  did  not 
cancel  out  as  it  happened  in  the  calculation  of  the 
equilibrium  tunneling  current  IQ.  A  density  of  states, 
rapidly  varying  with  energy,  could  then  produce  a 
structure  in  61.  However,  for  small  applied  bias, 
where  experimentally  the  blocking  current  is  observed, 
the  density  of  states  is  a  smooth  function  of  energy , 
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the  variation  of  which  would  not  be  expected  to  have 
significant  effect  on  61.  In  fact,  we  can  replace 


the  ratio 


D 


,  i  =  1,2,  by  its  value  at  the  Fermi 


27Th  Pj_ 

level  providing  eV  <<  ,  i  =  1,2 
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2iThp  . 


D(Efi> 

2L  . 
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fi 
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i=l ,  2 


(2.25) 


B  l 


where  V£ .  is  the  Fermi  velocity  and  L.  is  the  thickness 
f  l  i 

of  the  metal  i. 


t B  is  the  tunneling  lifetime  of  an 


.  th 


electron  in  i  electrode. 


Referencing  the  energies  to  the  Fermi  level  of 
the  left  hand  metal,  and  realizing  that  due  to  the 
difference  in  the  Fermi  functions  the  only  significant 
contribution  to  61  arises  from  energies  close  to  the 
Fermi  level,  we  rewrite  eg.  (2.24)  as 
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where 
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£/kBT 
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and  a  is  the  ideal  zero  bias  conductance  as  defined 
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in  eq.  (1.10) 
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Eq.  (2.26)  can  be  further  simplified  realizing  that  t, 

the  effective  relaxation  time  is  an  even  function  of 

the  excitation  energy  since  t  by  definition  (eq.  2.10) 

j  j.  ep 

is  even.  Using 


F°(e)  -  F°(e  +  eV)  =  -F°(-e-eV)  -F°(-s)  (2.28) 


we  rewrite  eq.  (2.26)  as: 
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The  blocking  conductance  6a (V)  =  6I(V)/dV,  normalized 
to  zero  bias  is  given  as: 


6a  (V) 
a 
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—  oo 


3F°  (s) 

3  £ 


r T  n  ( £  +  eV )  T  ^  ( £  +  eV )  ' 


(2.30) 


It  may  be  easily  verified  that  6a  (V)  is  an  even  function 
of  the  applied  bias  V.  Therefore,  it  is  symmetric  around 
zero  bias . 


In  further  discussions  we  assume  that  the  junc¬ 
tion  is  symmetric  and  that  the  applied  magnetic  field 
is  equal  to  0  kG .  We  write  then 


T  B  T  B1  1 B2 


and  t^(e)  =  T2^£)  “  T(£)  • 


At  0  K  we  can  evaluate  the  blocking  conductance  in  a 


closed  form, 
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Recalling  eqs .  (2 . 12 ) - (2 . 14 )  ,  we  write 
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for  a  F  = an<jjn /  n  =  1 ,  and 


6a  (V) 


2/t 


B 


a 


o 


2iTa 


+  (■ 


3fi 


2  s  |  ,7i  3 
— )  eV 


(2.34) 


for  n  =  2.  Figs. 2. 4  and  2.5  illustrate  the  blocking 
conductance  and  its  derivative  d6a (V) / dV  for  n  =  1  and 
n  =  2.  For  finite  temperatures  the  integral  in  eq . 
(2.30)  has  to  be  calculated  numerically.  We  write 
eq.  (2.30)  then  as 


6a (V)  = 


dx 


(ex  +  l)  (e  x+l) 


t ,  (xk„T  +  eV) 
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(2.35) 
d6a (V) 


The  derivative  of  the  blocking  conductance 
given  by 
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Blocking  conductance  and  its  derivative  at  T=0  K 

/  .2^  n  i  ^ 

(assuming  a  F  =  a^co  ,  n  =  1)  . 
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Fig .  2.5 


Blocking  conductance  and  its  derivative  at  T=0  K 
2 

(assuming  a  F  =  anw  ,  n  =  2) . 


39 


40 


1  d  6o(V)  _ 


a 


o 


dV 


(ex+l)  (e  x+l) 


[S, (xk_T  +  eV) 
_L  r> 


+  S2 (xkgT  +  eV) ] 


(2.36) 


where 
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To  parameterize  the  theoretical  curves  for  T  ^  0  K,  it 
is  convenient  to  introduce  a  parameter  r  ,  given  by: 
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j  =  ± 

n  ^x.aj 
i  n 


j  =  1,2 

n  =  1,2 


(2.38) 


This  follows  directly  from  the  limiting  behaviour  of  the 
effective  relaxation  time  at  low  temperatures. 
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Figs.  2.6  and  2.7  illustrate  the  behaviour  of 
the  blocking  conductance  at  zero  bias  for  a  symmetri¬ 


cal  junction  as  a  function  of  temperature.  A  more 
sensitive  comparison  of  the  theory  with  the  experi¬ 
ment  is  obtained  by  examining  the  behaviour  of  the 

1  d6a 


maxima  (minima)  in 


dV 


with  temperature.  Figs. 
1  da 


2.8  and  2.9  show  that  the  maximum  in  ^  decreases 
in  size  as  the  temperature  increases.  The  data 
calculations  were  normalized  to  1  K.  The  correspon¬ 
ding  peak  positions  shift  to  a  higher  bias. 

The  application  of  a  magnetic  field  leads  only 
to  a  renormalization  of  the  scattering  time  associated 
with  the  impurities. 
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2 

This  for  example,  at  T  =  0  K,  and  for  a  F  ^  m  ,  results 

in  a  shift  of  the  position,  V  of  the  maximum  in 

max 

1  d6a 


°o  dv  * 


Using  eg.  (2.33)  we  find 
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There  is  a  corresponding  change  in  the  peak  height: 
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Fig.  2.6 

Temperature  dependence  of  the  blocking 

2 

conductance  at  zero  bias  (a  F  =  a^oo)  . 


(arbitrary  units) 
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Fig .  2.7 


Temperature  dependence  of  the  blocking 

2  2 

conductance  at  zero  bias  (a  F  =  )  . 


(arbitrary  units) 
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Fig .  2.8 

Temperature  dependence  of  peak  positions, 

Vmax  (uPPer  diagram)  and  peak  amplitudes 
2_  d<5a 

- v  (lower  diagram)  for  different  values 

x>  av  max 

of  T  (a2F  =  a-j^oo )  . 
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Fig.  2.9 


Temperature  dependence  of  peak  positions, 


V  v  (upper  diagram)  and  peak  amplitudes 

lllclX 


_1 

a 


d6a 


dV 


V 


(lower  diagram)  for  different  values 


o  max 

of  F  (a  F  =  a.?a)  )  . 
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CHAPTER  3 


EXPERIMENTAL  TECHNIQUES 


3. 1  Sample  preparation 

The  preparation  of  tunnel  junctions  was  carried 
out  in  a  vacuum  system  by  an  evaporation  of  metals  onto 
a  suitable  substrate.  There  are  basically  three  steps 
involved  in  the  process  of  preparation.  First  a  metal 
film  (base  layer)  of  required  thickness  is  deposited. 

A  thin  insulating  layer  is  then  formed  on  the  top  of 
the  evaporated  film.  Finally  a  cross-strip  of  an  another 
metal  is  evaporated,  thus  forming  a  composite  metal- 
insulator-metal  (MIM)  system.  Evaporated  layers  of  In 
or  Ag  served  as  contacts  for  electrical  connection. 

The  junctions  were  prepared  on  Fisher  brand 
microscopic  glass  slides  cut  to  a  1/2  inch  square. 

The  slides  were  washed  with  a  detergent,  air  dried  and 
flame  polished.  For  some  Bi/Bi  junctions  freshly  cleared 
mica  was  used,  washed  in  methyl  alcohol.  The  substrates 
were  further  cleaned  by  exposing  them  to  the  plasma  of 
glow  discharge  (Holland,  1961)  prior  to  film  deposition. 
The  size  of  the  substrates  allowed  to  make  two  junctions 
with  a  common  base  layer  without  breaking  the  vacuum. 

The  films  were  evaporated  from  a  resistively 
heated  nonvolatile  support  (filament)  through  a  pattern 
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forming  mask,  in  close  vicinity  of  which  the  sub¬ 
strate  was  placed.  The  filaments  were  mounted 
on  a  rotating  carousel  which  allowed  a  proper 
indexing  of  the  evaporated  material.  Cylindrical 
quartz  chimneys  placed  around  them  prevented 
a  mutual  cross-contamination  of  the  metals.  Table 

3.1  lists  the  selection  of  filaments  for  evaporation. 

The  distance  between  the  vapour  source  and  the  sub¬ 
strate  was  7  inches. 

The  vacuum  chamber,  enclosed  in  an  18  inch  bell 
jar,  was  pumped  by  an  oil  diffusion  pump  with  a  liquid 
nitrogen  cold  trap  backed  by  a  mechanical  pump.  A 
zeolite  trap  was  installed  between  the  mechanical  pump 
and  the  rest  of  the  system.  The  system  was  ready  for  the 
evaporation  after  1  hour  of  evacuation  with  a  pressure 
lower  than  2x10  ^  torr . 

3 . 2  Formation  of  tunnel  barriers 

A  glow  discharge  technique  (Miles  and  Smith,  1963) 
was  used  to  prepare  tunnel  junctions  of  Mg  and  A1 .  It 
consists  of  setting  up  a  discharge  between  two  aluminum 
electrodes  while  the  pressure  of  the  system  was  maintained 
between  100-600  torr  of  either  99.5%  pure  oxygen  or 
nitrogen  of  the  same  purity  (for  Mg  only) .  The  aluminum 
electrodes  were  2  inches  wide,  separated  1/2  inch 


- 


' 


54 


Table  3.1 

Selection  of  Materials  for  Evaporation 


Metal 

Purity 

Supplier 

Evaporation  Source 

Al 

59 

A . D .  Mackay 

F8-3X. 030  W 

coil 

Mg 

59 

Johnson  Mattey 

S17A-.005  Ta 

covered  boat 

Bi 

69 

Cominco 

S39-.005  Ta 

S39-.005  Mo 

boat 

Pb 

69 

Cominco 

S39-.005  Mo 

boat 

Au 

59 

Cominco 

S39-.005  Mo 

S9A-.005  W 

boat 

Evaporation  sources  supplied  by  R.D.  Mathis  Company 
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apart  and  placed  2  inch  above  the  holder  containing 
the  glass  slide  with  a  freshly  evaporated  metal  film. 
The  polarity  of  the  applied  DC  voltage  was  such  that 
the  evaporator  which  was  conductively  connected  with 
one  of  the  electrodes  was  at  positive  potential  with 
respect  to  the  second  electrode.  The  thickness  of  a 
tunnel  barrier  grown  by  this  method  depended  on  the 
time  the  discharge  was  maintained,  pressure  in  the 
system  and  the  current  flowing  through  the  electrodes. 
For  a  successful  preparation  of  a  tunnel  barrier  it  was 
important  to  keep  a  stable  discharge,  without  random 
fluctuations.  Presumably  the  rapidly  varying  electric 
field  produced  pinholes  in  the  thin  insulator.  The 
oxidation  of  Mg  for  1-3  minutes  in  100-200  torr  of  0^, 
with  a  glow  discharge  current  of  5-15  mA  was  sufficient 
to  produce  tunnel  junctions  with  resistances  of  nomin¬ 
ally  200  fi.  A1  required  a  longer  oxidation  time 
(5-10  minutes)  and  higher  pressure  (400  torr) ,  while 
the  current  was  5-15  mA.  To  prepare  junctions  in  ni¬ 
trogen  environment,  the  discharge  had  to  be  kept  for 
40-80  minutes  with  currents  up  to  40  mA  and  pressures 
ranging  between  400-600  torr. 

Repeated  experiments  of  producing  an  oxide  of 
Bi  by  the  glow  discharge  technique  ended  in  failure. 
However,  a  method  of  thermal  oxidation  (Hapase  et  al, 
1967)  yielded  thick  enough  oxides  to  produce  junctions 
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having  resistances  around  200  ft.  The  oxidation  procedure 
was  as  follows:  After  the  deposition  of  the  Bi  film 
the  vacuum  was  broken  and  the  substrate  was  placed  in 
an  oven  at  180°-220°C  for  2  hours.  Following  this,  the 
substrate  was  again  placed  in  the  substrate  holder  and 
the  vacuum  chamber  re-evacuated  for  the  evaporation  of 
the  cross-film. 


3 . 3  Film  thickness  determination 

The  thickness  of  films  during  vacuum  depo¬ 
sition  was  monitored  by  Sloan  DTM-3  (Sloan  Instr.  Corp.) 
quartz  crystal  monitor.  The  measurement  utilizes  a 
principle  that  the  resonance  frequency  of  a  piezoelectric 
crystal  depends  on  the  mass  of  the  material  deposited 
on  the  surfaces  of  the  crystal  (Eschbach  and  Kruidhoff, 
1965).  The  monitor  therefore  consists  of  a  crystal 
placed  adjacent  to  a  substrate  onto  which  the  film  is 
deposited  and  the  electrical  equipment  for  measuring 
the  change  in  its  resonance  frequency.  Eq .  (3.1) 

relates  the  thickness  L  of  the  deposited  metal  to  the 
measured  frequency  shift  Af  ( Stockbridge ,  1966) 


L  =  C1C2 


Af 


CAf 

P 


(3.1) 


C  =  C1C2 


is  a  constant  dependent  on  the  thickness  and  cut  of 
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the  piezoelectric  crystal,  C2  depends  on  the  geometry 
of  the  distribution  of  the  deposit  over  the  crystal 
and  p  is  the  density  of  the  material.  The  validity  of 
eq .  (3.1)  is  restricted  to  mass  deposits  which  are  at 

most  5%  of  the  mass  of  the  crystal  itself  (Steckelmacher , 
1965) .  For  all  practical  purposes  this  condition  was 
satisfied.  Since  the  density  of  thin  films  is  different  from 
that  of  the  bulk  due  to  voids  in  the  film  (Hartman,  1965) 
it  is  best  to  determine  the  mass  sensitivity  ratio  C/p 
for  each  given  material  experimentally,  combining  the 
shift  in  the  resonance  frequency  of  the  crystal  with  an 
optical  thickness  measurement  using  Tolansky  interferro- 
meter  (Tolansky,  1960). 

The  interferometer  consists  of  two  slightly 
inclined  optical  flats,  one  of  them  supporting  the  film 
which  forms  a  step  on  the  substrate.  When  the  second 
optical  flat  is  brought  in  contact  with  the  film  surface 
and  the  interferometer  illuminated  with  a  parallel 
monochromatic  beam  at  normal  incidence,  dark  fringes 
are  observed  which  trace  out  the  points  of  equal  air- 
gap  thickness.  Their  separation  is  A/2  where  A  is  the 
wavelength  of  the  monochromatic  light.  The  fringes 
show  a  displacement  as  they  pass  over  the  film  step 
edge.  The  displacement  d  expressed  as  a  function 

of  A/2  fringe  spacing  gives  the  film  thickness. 

_ d _  A 

fringe  spacing  2 


L 


(3.2) 
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In  order  that  the  phase  change  on  reflection  from  the 
two  sides  of  the  step  will  be  the  same,  it  is  necessary 
to  coat  the  film  as  well  as  the  exposed  glass  surface 
with  the  same  reflecting  silver  layer. 

Vrba  (Vrba,  1971)  has  described  these  measure¬ 
ments  in  greater  detail.  The  experimentally  determined 
values  of  C/p  are  the  following: 


Al  - 

940 

A/k 

Mg  - 

1480 

A/k 

Bi  - 

268 

A/k 

Pb  - 

235 

A/k 

Au  - 

127 

A/k 

The  uncertainties  in  the  mass  sensitivity  determination 
are  approximately  5%  and  are  of  interferometric  and 
geometric  origin. 

3 . 4  Temperature  measurement  and  control 

The  temperatures  were  determined  by  a  combination 
of  vapour  pressure  and  germanium  resistance  thermometer 
measurements.  A  1958  scale  of  temperatures  (Dijk  et  al , 
1960)  was  used  to  convert  the  vapour  pressure  to  tem¬ 
perature.  Above  the  boiling  of  He^ ,  a  germanium  ther¬ 


mometer  (Scientific  Instruments,  #  10372)  was  used  with 
a  computer  generated  polynomial  fit  of  the  form 
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b  =  l  akHk(ln  R)  (3.3) 

T  k=0  K  K 

where  ak  is  the  kth  expansion  coefficient,  Hk  -  the 
Chebyshev  polynomial  of  order  k,  T—  temperature,  R. 
resistance  of  the  thermometer. 

4 

By  a  controlled  regulation  of  the  He  vapour 
pressure,  temperatures  between  1  K  and  4.2  K  were 
obtained.  Sommers  (Sommers,  1954)  described  a  con¬ 
venient  manostat  which  employs  a  thin  teflon  diaphragm, 
which  bows  out  when  there  is  a  pressure  difference 
across  it,  thus  blocking  the  pumping  orifice.  The 
temperature  fluctuations  using  this  manostat  did  not 
exceed  5  mK  around  a  preset  value.  A  solid  state 
temperature  controller  designed  by  Dr .  Adler  and  Dr . 
Rogers  was  used  to  regulate  temperatures  above 
4.2  K.  The  circuit  is  of  the  conventional  Wheat¬ 
stone  bridge  type,  using  a  germanium  thermometer  as  its 
sensing  element  and  providing  a  proportional  heater 
feedback  to  the  regulated  bath.  It  was  capable  of 
maintaining  constant  temperature  within  ±1  mK  at  5  K. 
The  power  dissipation  in  the  thermometer  was  less  than 
10-8  W. 
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Fig .  3.2 


Schematic  diagram  of  the  cryostat 
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3 . 5  Production  of  low  temperatures 

4 

A  conventional  He  cryostat  facilitated  measure¬ 
ments  on  tunnel  junctions  at  low  temperatures.  A 
schematic  diagram  of  the  cryostat  is  illustrated  on 
fig.  3.2.  Indium  0-rings  with  Apiezon  N  grease  were 
used  to  seal  the  vacuum  and  sample  chambers.  The 
cryostat  was  designed  to  fit  into  the  bore  of  a  70  kG 
superconducting  magnet.  By  pumping  on  the  helium  liquid 
level,  in  the  sample  chamber,  temperatures  as  low  as 
1  K  were  attained. 

3 . 6  Measurement  of  tunneling  characteristics 

A  system  using  a  combination  of  harmonic  detection 
and  bridge  techniques  (Adler  and  Jackson,  1966) 
suitably  interfaced  to  an  online  mini-computer  facility 
for  an  analog  and  digital  recording  was  used  in  obtaining 
the  relevant  information  on  tunnel  junctions.  These 
are  the  dynamical  conductance  a  =  dl/dV  and  its  deriva- 
tive  da/dV  =  d  I/dV  .  The  principle  of  harmonic  detec¬ 
tion  method  can  be  understood  by  considering  fig.  3.3. 

In  this  circuit  the  tunnel  junction  serves  as 
one  arm  of  an  ac  Wheatstone  bridge,  the  other  arm  be¬ 
ing  a  variable  resistor  Rd.  Assuming  that  the  bridge 
is  close  to  its  balance  the  current  through  its  arms  is 
nearly  identical.  In  the  voltage  region  where  the 
tunnel  junctions  are  investigated,  they  behave  as  weakly 
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Fig .  3.3 


Simplified  bridge  configuration 
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nonlinear  passive  elements.  If  the  small  current 
modulation  amplitude  6  is  kept  constant,  then  the 
voltage  V  across  the  junction  can  be  written  in  terms 
of  a  Taylor's  series 

2 

V  ( I )  =  V  ( I  )  +  (™)  6  cos  w  t  +  ~  (— y)  S2cos2cot 

°  dI  In  2  dl  i 

O  -Lo 

+  .  (3.4) 


or 

V(I) 


V(Io) 


(dV) 

vdl;  x 


6  cos  cot  + 


1 

4 


o 


<^> 

dl 


(1+2  cos  2cot) 


'O 


+ 


(3.5) 


where  V(I  )  is  the  DC-bias  across  the  junction  and  go  is 
the  modulating  angular  frequency.  The  voltage  across 
the  variable  resistor  Rd  which  has  a  linear  I-V 
characteristics  is  given  by 


Vd  =  (I  +  6  cos  oot)  Rd 


(3.6) 


Thus  the  difference  in  potential  between  points  1  and 
2 ,  V^2  r  has  components  at  the  modulating  frequency  m  and 
its  various  harmonics.  The  signal  at  the  fundamental  fre¬ 
quency^^  (w),  is  proportional  only  to  the  departure  of  the 
dynamic  resistance  (^0  t  f rom  an(3  the  second  harmonic 


d2V 


Both 


component,  V-^  (2co)  ,  is  proportional  to  * 

these  components  are  measured  by  means  of  a  phase  sensitive 
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detector . 


V12(o))  =  <5[Rd  -  (gj)  ]  cos  cot 

^o 


(3.7) 


V19(2co)  =  j  6^  (  — y)  cos  2 cot  .  (3.8) 

dl  I 

o 

A  direct  calibration  of  curve  enables  the 

dl 

computation  of  the  dynamical  conductance  a  = 


Ideally,  the  second  derivative  — i  ^ 

dv2  dv 

culated  using  the  identity: 


is  then  cal- 


2  2 
d  I  3  d  V 

- =  -a  •  - ~- 


(3.9) 


dv' 


dl' 


The  presence  of  a  second  harmonic  distortion  in 
the  modulation  supply  complicates  the  above  relation. 
This  necessitates  a  use  of  a  different  numerical  method 
in  calibrating  The  mathematical  details  of  the 

analysis  of  this  problem  are  given  in  Appendix  I  with 
an  earlier  version  of  our  data  acquisition  system.  The 
present  system  which  employs  an  online  mini-computer  and 
its  peripherals  allows  a  fast  evaluation  of  the  incoming 
data.  The  complete  description  of  the  new  data  acquisi¬ 
tion  system  was  given  in  a  technical  report  (Adler, 
Straus ,  1974). 


. 


; 


Fig .  3.4 


Schematic  diagram  of  the  data  acquisition 


system 


CHAPTER  4 


EXPERIMENTAL  RESULTS 

4 . 1  Introduction 

To  test  the  theory  with  the  experiment  we  have 
prepared  a  series  of  tunnel  junctions  using  different 
metals  as  their  electrodes.  The  thickness  of  these 
ranged  between  500-50000  A.  The  insulating  layers 
of  the  junctions  were  prepared  either  by  oxidation 
or  nitrification  (for  Mg  only) .  Some  junctions 
were  deliberately  doped  with  oil  to  reduce  the  pro¬ 
bability  of  making  a  sample  with  spurious  oscillations 
in  their  characteristics.  Junctions  having  these,  will 
be  discussed  in  a  separate  chapter. 

4 . 2  The  Problem  of  the  Background 

A  typical  measurement  of  the  conductance  of  a 
Mg/oil/Mg  junction  to  a  5  meV  range  is  shown  in  fig. 

4.1  for  several  different  temperatures.  As  the  tem¬ 
perature  increases,  the  conductance  minimum  decreases 
with  a  simultaneous  loss  in  the  sharpness  of  the 
anomaly.  The  second  derivative  measurement  for  the 
same  junction  is  shown  in  the  next  figure.  Here, 
the  maximum  of  da/dV  with  higher  temperatures  decreases 
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Fig.  4.1 


Temperature  dependent  conductance 
for  a  Mg/Mg  junction. 
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Fig.  4.2 


Temperature  dependent  normalized 
derivative  of  the  conductance  for 
a  Mg /Mg  junction. 
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(for  a  positive  applied  bias)  while  its  position  shifts 
to  a  higher  voltage. 

The  conductance  curve  in  fig.  4.1  has  a  slight 
asymmetry  about  the  origin.  We  assume  that  this 
asymmetry  arises  from  the  ideal  background  conductance 
as  discussed  earlier  (section  2.1).  To  eliminate  the 
background,  which  we  consider  to  be  temperature 
independent  at  least  between  1-10  K,  we  have  adopted 
the  following  two  procedures. 

Let  us  first  recall  fig.  2.6,  where  we  notice 
that  at  higher  temperatures,  the  dip  in  conductance 
due  to  blocking  is  essentially  equal  to  zero.  There¬ 
fore  by  subtracting  the  conductance  measured  at  the 
highest  temperature  (10  I<)  from  the  ones  measured  at 
lower  temperatures  where  the  blocking  is  still  effec¬ 
tive,  we  can  isolate  the  pure  "blocking  conductance 
6a".  This  technique  was  successfully  applied  for  the 
junction  shown  earlier.  Fig.  4.3  and  fig.  4.4  illus¬ 
trate  the  resulting  curves. 

For  junctions  where  the  high  temperature  curve 
was  not  obtained  we  have  considered  the  following 
alternative:  Let  us  suppose  that  the  experimental 

conductance  may  be  represented  by  an  expression  of 
the  form 

a (V)  =  a  +  a.V  +  a„V2  +  Sa(V,T) 
o  l  z 


(4.1) 
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Fig .  4.3 


Temperature  dependence  of  the  blocking 
conductance  for  a  Mg/Mg  junction  after 
subtracting  the  background. 
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Fig.  4.4 


Temperature  dependence  of  the 
normalized  derivative  of  the  blocking 
conductance  for  a  Mg/Mg  junction  after 
subtracting  the  background. 
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where  the  first  three  terms  constitute  the  ideal 
temperature  independent  background.  Then  by  taking 
the  even  part  of  the  conductance 


0e(V) 


O 

O 


oo 


(g(V)  +  g(-V)  _  2  +  6_g 

v  2  o  g 

o  o 


(V]_T)  +  1  . 


(4.2) 


We  eliminate  the  linear  term.  The  assumption  that  g-^ 
is  temperature  independent  is  verified  by  forming 
the  odd  part  of  the  conductance 


n  =  Q  (V)  -  0  (~v) 
odd  2 


°iv 


(4.3) 


and  plotting  it  as  function  of  the  bias  voltage  for 

different  temperatures.  The  resulting  curve  should 

be  a  straight  line  independent  of  the  temperature. 

This  really  is  the  case,  as  fig.  4.5  demonstrates. 

2 

To  eliminate  the  quadratic  term  V  ,  a  more  detailed 
knowledge  about  the  background  is  needed  which  unfor¬ 
tunately  we  do  not  have.  We  can,  however,  use  the 
fact  that  the  blocking  at  4-5  meV  is  negligible  because 
of  the  short  lifetimes  of  the  excitations.  Therefore 
it  is  appropriate  to  subtract  from  the  even  part  of  the 
second  derivative 


dg 


e 


d6g  (V ,T ) 


dV 


2g2V  + 


dV 


(4.5) 
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Fig .  4.5 


Experimental  odd  conductance 


odd 
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a  line ,  passing  through  the  origin,  with  a  slope  which 
brings  the  measured  da^/dV  to  zero  at  4-5  meV.  We 
estimate  that  this  technique  may  introduce  at  the  most 
a  5%  uncertainty  in  the  determination  of  the  peak 
magnitudes  in  dae/dV. 

4.3  General  Considerations  in  Extracting  x^,  x^  and 
Tdrif t 

A  successful  theory  of  the  zero  bias  anomaly 
should  be  able  to  predict  the  lineshape  of  the  conduc¬ 
tance  dip  and  its  general  behaviour  with  temperature 
or  magnetic  field. 

Let  us  consider  first  the  case  of  a  zero  magni- 
tic  field.  The  usual  procedure  in  comparing  the 
experimental  data  with  the  theory  is  to  fit  the  mea¬ 
sured  curve  at  a  fixed  temperature  and  extract  the 
necessary  parameters  characterizing  the  anomaly.  These 
are  x^  -  the  tunneling  time  and  x^  -  the  impurity 

scattering  time.  We  assume  that  the  coupling  strength 
of  the  electron  phonon  interaction  is  known  beforehand. 
Having  obtained  x^  and  x^,  the  theoretical  curves  are 
then  calculated  for  different  temperatures  and  compared 
to  the  measured  ones.  In  practice,  it  is  more  sensi¬ 
tive  and  accurate  to  fit  the  zero  bias  anomaly  to  the 
experimentally  determined  derivative  of  the  even  part 
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,  d6a 

of  the  conductance  —  — The  curve  is  fitted  with 

o 

respect  to  the  peak  height  and  its  peak  position  in 

terms  of  Y  and  the  ratio  of  T^/Tg*  For  t^ie  sa^e  oF 

completeness  we  rewrite  eq.  (2.30)  in  these  terms 

2 

(assuming  that  a  F  =  a^u)) 

CO  f 


1 


a 

o 


2x . .  r 


_ dx 

(ex+l)  (e"X+l) 


_ 1 _ 

k  (F1+  (xkgT+eV)  2+  (TTkBT)2) 


+ 


1 


(f1+  (eV-xkfiT) 2  +  (rkfiT)2J 


(4.6) 


The  derivative  of  even  conductance  is  given  as 


d6a  4t  .  r, 
_1 _ e  _  l 

a  dv 
o 


B 


dx 


0 


(ex+l)  (e  x+l) 


(xk„T+eV) 

i D 


(T1+  (xkgT+eV ) 2+ (rkgT) 2) 2 


+ 


(eV-xkBT) 


(1^+  (eV-xkgT )  2+  (rkBT)  2)  2J 


(4.7) 


When  writing  the  above  equations  we  have  again 
assumed  that  the  junction  is  symmetric  with  respect  to 
its  relaxation  processes  and  tunneling  times.  ihis 
simplifying  assumption f as  well, was  used  in  the  analysis 
of  the  experimental  data.  From  the  theoretical  point 
of  view ,  it  is  hard  to  justify  this,  even  for  samples 
prepared  with  the  same  electrodes,  having  the  same 
Generally,  the  films  can  have  different 


thickness . 


. 


34 


grain  sizes  and  imperfections.  In  addition,  the  two 
interface  planes  adjoining  the  insulating  layer  can 
be  contaminated  with  different  impurities.  This  all, 
leads  to  a  tunnel  junction  with  non- symmetrical  tunnel¬ 
ing  and  scattering  times.  Considering  eq.  (2.30)  we 
should  then  write  the  part  enclosed  in  brackets  as 


Bl 


(r];+£2+(7TkBT)2) 


+ 


Bl 

B2 


1 

2 


al  (ri+e2+(7TkBT)2) 


(4.8) 


from  which  it  follows  that  generally  we  can  use  three 
1  2 

parameters  r^,  and  A 


T 


A  = 


Bl 

B2 


a 

a 


1 

1 

2 

1 


(4.9) 


to  fit  the  temperature  dependence  of  the  experimental 
data.  The  fourth  parameter  xB1  adjusts  the  peak  ampli¬ 
tudes.  With  a  three  parameter  fit  we  have  to  consider 
the  interference  between  the  two  terms  in  (4.8)  which 

can  be  responsible  in  creating  an  additional  peak  in 

.  1  d6ae 

dSae/dv.  Experimentally,  the  two  peaks  m  dv  were 

never  observed.  This  does  not,  however,  exclude  the 
possibility  of  their  insufficient  separation  and  the 

small  size  of  the  second  peak. 

Numerical  analysis  (Teshima,  1974)  indicates 
that  the  second  derivative  will  not  have  the  two  peaks 


as  long  as 
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<  8  +  /90 


at  T  =  0  K  . 


Since  the  experimentally  found  parameter  F  ranged  from 

2 

.4  <  T  <2. (meV)  ,  it  is  safe  to  assume  that  the 
condition  above  is  always  satisfied. 

Nevertheless  we  have  compared  some  of  our  data 
with  the  predictions  of  the  theory  using  a  three  para¬ 
meter  fit.  In  general,  there  was  no  appreciable 
difference  between  this  fit  and  the  one  obtained  by 
assuming  a  symmetrical  junction. 

Furthermore,  we  assume  that  t  ,  the  relaxation 

ep ' 

time  due  to  electron-phonon  interaction  is  given  by 

eq.  (2.12).  This  is  equivalent  to  an  assumption  that 
2 

a  F  behaves  linearly  for  small  excitation  energies. 

The  justification  of  this  assumption  is  partly  based 
on  the  preliminary  results  of  measurements  on  Mg/Mg 
and  Mg/Au  junctions  down  to  0.3  K.  The  second  deriva¬ 
tive  characteristics  of  these  junctions  have  not  shown 

kink  in  the  vicinity  of  zero,  whereas  they  would  if 
2 

a  F  behaved  quadratically  for  small  energies  (fig.  2.5). 

In  addition,  the  theoretical  fits  using  t  derived 

^  ep 

2 

on  the  basis  of  the  linear  a  F  behaviour  were  superior 

to  those  obtained  by  t  of  eq .  (2.13)  and  eq.  (2.14). 

ep 

The  relaxation  time  x,  . - ,  associated  with  the 

drift 

magnetic  field  measurements  can  be  extracted  only  if 
both  x.  and  x„  are  known.  But  these  are  obtained  from 
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zero  field  measurements.  The  procedure  is  strictly 
valid  for  T-  which  is  magnetic  field  independent.  The 

ID 

results  prove  this  assumption  to  be  a  correct  one. 

4.4.  Comparison  of  the  Theory  with  the  Experiment 

Let  us  first  consider  the  results  of  the  theory 

applied  on  the  previously  discussed  Mg/Mg  junction. 

in  fig.  4.6,  beside  the  experimentally  determined 

data  points,  we  show  the  behaviour  of  the  theoretical 

curves  for  three  different  parameters  T  =  0.4,  0.8, 

1.2  meV  .  The  peak  amplitudes  were  normalized  at  1  K. 

The  dependence  of  the  peak  positions  on  the  temperature 

2 

is  well  fitted  by  a  theoretical  curve  with  F  =0,4  meV  . 

In  this  case,  however,  the  peak  amplitudes  do  not  drop 

as  fast  as  the  theory  would  predict.  The  second  extreme, 

T  =1.2  meV2,  fits  the  peak  amplitudes  better  at  higher 

temperatures,  whereas  the  peak  positions  lie  below  the 

theoretical  curve.  A  reasonable  compromise  is  to  assume 

2 

that  r  =  0.8  meV  represents  the  best  fit. 

The  theory  then  predicts  a  stronger  temperature 
dependence  of  the  peak  amplitudes  than  experimentally 
observed.  Correspondingly  the  peak  positions  are  at 
lower  biases  than  expected. 

2 

The  reasons  behind  accepting  T  =  0.8  meV  as  the 
best  fit  will  be  more  transparent  when  we  consider  the 
results  of  the  experiments  on  other  junctions.  We 
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Fig.  4.6 


Comparison  of  the  experimentally  determined 


peak  positions ,  (upper  diagram)  and  peak 

max 

d6oQ 

■  (lower  diagram) 

max 


amplitudes  , 


1 


dV 


for  a  Mg/Mg  junction  at  various  temperatures 
with  peak  positions  and  peak  amplitudes 
predicted  by  the  theory. 
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will  find  that  these  conclusions  apply  equally  for 
them . 

2 

With  T  =0.8  meV  we  have  plotted  the  tempera¬ 
ture  dependence  of  the  blocking  conductance  at  zero 
bias  and  compared  it  with  the  experiment  (fig.  4.7) 
the  data  were  normalized  at  1  K. 

To  extract  t.  and  xD  we  have  to  know  the  value 
of  the  electron-phonon  coupling  constant  for  Mg. 

Since  to  our  knowledge  there  is  no  available  informa¬ 
tion  on  it,  we  have  taken  it  to  be  equal  to  .0015  meV 
This  value  is  one  half  of  that  published  for  aluminum 
(Carbotte,  Dynes,  1968).  This  certainly  introduces 
an  uncertainty  in  the  determination  of  t.  and  td. 
Nevertheless,  we  feel  that  this  choice  is  appropriate 
in  obtaining  the  order  of  magnitudes  involved. 

We  have  then 
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Using  these  parameters  we  can  compare  the  theoretical 
even  part  of  the  second  derivative  with  the  measured 
one.  Fig.  4.8  illustrates  the  fit  at  T  =  1  K  and 
3  K.  The  peak  amplitude  at  3  K  reproduces  well,  while 
the  experimental  peak  position  is  shifted  to  a  lower 


value . 
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Fig.  4.7 

Comparison  of  the  temperature  dependent 
blocking  conductance  at  zero  bias  for  a  Mg/Mg 
junction  with  the  blocking  conductance  cal¬ 
culated  from  the  theory. 
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Replacing  the  top  Mg  layer  by  an  Au  counter¬ 
electrode  has  not  altered  the  character  of  the  zero 
bias  anomaly.  Figs.  4.9  and  4.10  show  the  temperature 
dependence  of  the  conductance  of  two  Mg/N^/Au  junctions. 
These  were  prepared  on  the  same  substrate  with  the  same 
Au  overlayer.  Completing  their  initial  measurement, 
the  junctions  were  allowed  to  warm  up  while  they  still 
were  kept  under  a  vacuum.  After  a  period  of  a  few  days, 
the  junctions  were  again  cooled  down  and  remeasured.  The 
annealing  produced  a  2.5  fold  increase  in  the  resistance 
of  the  tunnel  junctions.  The  relative  magnitude  of  the 
conductance  dip,  characteristic  of  the  zero  bias  anomaly, 
however,  has  not  changed.  This  documents  that  the 
relevant  scattering  mechanisms  which  deblock  the 
tunneling  electrons  were  not  altered.  In  fig.  4.11 

we  show  the  experimental  data  for  85  B1  and  B2  (annealed) 

2 

together  with  the  fit  obtained  by  adjusting  T  to  0.6  meV  . 
The  theory  again  predicts  a  stronger  temperature  depen¬ 
dence  for  the  peak  amplitudes.  As  the  temperature 

decreases  the  peak  positions  lie  below  the  theoretical 

1  d6ae 

curve.  The  lineshapes  of  — —  for  T  =  1  K  and  3  K 

°o  uv 

are  shown  in  the  following  figure. 

According  to  the  theory,  metals  with  infinitely 
thick  electrodes  are  essentially  in  equilibrium.  This 
is  because  the  density  of  final  states  for  the  incoming 
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Fig .  4.9 


Temperature  dependence  of  the  zero  bias 
conductance  anomaly  for  085A  Mg/N2/Au 

Mg  =  30000  A 
Au  =  2200  A 
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Fig.  4.10 

Temperature  dependence  anomaly  for  085B 
Mg/N2/Au 

Mg  =  3000  A 
Au  =  2200  A 
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Fig.  4.11 


Temperature  dependence  of  peak  positions,  vmax 

1  d(5ae 

(upper  diagram)  and  peak  amplitudes ,  — 

o  vmax 

(lower  diagram)  of  085,  before  (Bl)  and  after 
(B2)  annealing.  The  solid  lines  represent  the 
theoretical  fit. 
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Fig .  4.12 


Comparison  of  experimental  and  theoretical 
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electrons  is  infinite.  Mathematically,  this  is 
expressed  by  setting  the  tunneling  time  x  to  infinity. 

To  test  this  assertion  we  have  prepared  double 
junctions  which  had  one  common  electrode  and  in  which 
the  counterelectrodes  of  the  same  material  had  different 
thicknesses.  Samples  so  prepared  consistently  had 
smaller  zero  bias  anomalies  in  junctions  with  the 
thicker  film.  Junctions  85A,  85B  are  the  representative 
of  these.  The  Au  film  was  2200  A  thick.  The  Mg  elec¬ 
trode  for  the  junction  85B  was  3000  &  thick  while  for  85A 
it  was  30000  A  thick  (fig. 4. 13).  The  solid  lines  represent 

the  theoretical  fits.  The" thinner" junction  85B  was 

2  .  . 

parameterized  with  F  =  0.6  meV  ,  giving 
t„  =  1.3  X  10  8  sec 

ID 

x .  =  2.0  x  10  ^  sec . 

l 

2 

The  "thicker"  junction  was  fitted  w7ith  F  =  1.3  meV  yield¬ 
ing 

td  =  1.0  x  10  8  sec 

i D 

=  1.1  x  10  10  sec. 

Of  course,  these  tests  can  give  only  qualitative  answers. 
Meaningful  quantitative  numbers  could  be  obtained  only 
for  metal  electrodes  with  a  character  of  a  single  grain 
throughout  the  entire  film  thickness.  In  practice,  the 
evaporated  films  may  grow  in  layers  of  different  thick- 


. 


Fig .  4.13 


Zero  bias  blocking  conductance  for  two  junctions 
with  different  thickness  of  metal  electrodes. 

( 8 5 A  -  Mg  =  30000  A,  Au  =  2200  A) 

(85B  -  Mg  =  3000  A ,  Au  =  2200  A) . 

The  theoretical  fit  (solid  line) for  085B  was 
calculated  with  T  =  0.6  meV2 .  85A  was  fitted 

with  1  =  1.3  meV^ . 
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ness  and  orientation.  This  then  reduces  the  overall, 
observable  effect.  The  results  of  this  discussion 
apply  equally  well  to  that  part  of  the  impurity  rela¬ 
xation  time  tv  which  is  responsible  for  the  diffuse 
boundary  scattering.  An  investigation  of  Mg  films  by 
an  electron  microscope  revealed  that  Mg  crystallizes 
in  grains,  the  thickness  of  which  could  be  equal  to 
the  thickness  of  the  whole  specimen.  This  was  probably 
true  of  films  less  than  1500  A  (Turner,  1974) .  The 
lateral  dimensions  of  the  grains  ranged  between  500-800 

O 

A.  The  decrease  of  the  zero  bias  anomaly  in  junctions 
with  thicker  (>1500  A)  films  then  shows  that  a  fraction 
of  the  incoming  electrons  feels  the  increased  number  of 
grains  grown  across  the  whole  film. 

Similar  results  have  been  obtained  on  Mg/Au 
junctions  having  an  oxide  as  its  insulating  layer. 

In  fig.  4.14  we  show  a  fit  for  such  a  junction. 

The  zero  bias  anomaly  measurements  on  Mg/Au  and 
Mg/Mg  systems  have  been  complemented  with  investigations 
on  junctions  having  different  electrodes  (Pb,  Al) .  We 
illustrate  the  results  on  Al/Al  structures.  Fig.  4.15 
shows  a  temperature  dependence  of  the  conductance  dip  out 
on  a  previously  published  Al/Al  junction  (Trof imenkof f , 
et  al,  1972).  The  data  were  fitted  to  the  3  K  curve  width 
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Fig.  4.14 


Comparison  of  the  experimentally  determined  peak 
positions,  V  (upper  diagram)  and  peak  amplitudes, 


d6a 

_1 _ e 

a  dV 
o 


(lower  diagram)  for  a  Mg/O/Au  junction 


V 


max 


at  various  temperatures  with  peak  positions  and 
peak  amplitudes  predicted  by  the  theory. 
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Fig 


4.15 


Temperature  dependent  conductance 
of  an  Al/Al  junction. 
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(Carbotte,  Dynes,  1968)  .  Figs.  4.16,  4.17  show  the 
resulting  fit.  In  these  measurements  care  had  to  be 
taken  to  separate  the  decrease  in  conductivity  due  to 
the  blocking  effect  from  the  one  arising  from 
the  appearance  of  superconducting  energy  gap  in 
the  electrode.  This  is  particularly  true  of  thin 
films,  where  it  is  known  that  the  transition  tem¬ 
perature  increases  with  decreasing  film  thickness. 

We  show  this  effect  in  fig.  4.18,  where  at  2  K  the 
Al  film  (680  A  thick)  is  already  superconducting. 

The  transition  temperature  of  the  bulk  Al  is  1.19  K. 

In  principle  a  magnetic  field  could  have  been  used 
to  suppress  the  appearance  of  the  superconducting 
phase.  But  since  we  wanted  to  eliminate  the 
influence  of  the  magnetic  field  on  the  size  of  the 
zero  bias  anomaly  for  the  subsequent  field  dependence 
measurements  this  has  not  been  done. 

The  results  in  the  presence  of  the  magnetic 
field  agreed  with  the  predictions  of  the  theory 
(Sec.  2.5).  They  have  been  described  earlier  (Adler, 
et  al,  1973)  and  are  attached  as  an  Appendix  (Appendix 


II)  . 


i 


- 


112 


Fig.  4.16 


Peak  height  (open  circles)  and  peak  position 
(full  circles)  vs  temperature  for  an  Al/Al 
junction.  The  solid  line  represents  the 


theoretical  fit. 
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Fig.  4.17 


da 

e  TT 

dv~vs  v 

represent 


for  an  Al/Al  junction.  Solid  lines 
data;  dashed  lines  theoretical  result 


for  T 


3  K. 
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Fig.  4.18 


Illustration  of  the  appearance  of  the 
superconducting  energy  (at  2  K)  gap  m 
an  Al/Al  junction. 


851  AIAI 


117 


<N 


[on UJ  ]  D  xr0l 


CN 


ENERGY  meV 


118 


Fig.  4.19 


Temperature  dependence  of  the  blocking 
conductance  at  zero  bias  for  691  Al/Al  and 
851  Al/Al.  The  solid  lines  are  just  indicating 


the  linear  behaviour. 
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4 . 5  Discussion 

In  the  last  section  we  have  presented  the  results 
of  the  experimental  investigations  on  small  conductance 
decreases  around  zero  bias.  These  were  compared  with 
the  predictions  of  the  transport  model,  based  on  the 
non-equilibrium  tunneling  theory.  The  temperature 
dependence  of  the  anomaly  was  fitted  by  one  parameter, 
the  knowledge  of  which  allowed  to  extract  the  impurity 
relaxation  time,  xi,  if  the  strength  of  the  electron- 
phonon  coupling  was  known  beforehand.  The  amplitude  of 
the  theoretical  curve  was  adjusted  to  the  experimental 
one  with  an  aid  of  another  parameter,  tb,  related  to  the 
barrier  penetrability  coefficient  and  the  electrode 
thickness . 

The  theory,  as  it  stands,  predicts  a  somewhat 
stronger  temperature  dependence  of  the  observed  peak 
amplitudes  in  dae/dV.  The  peak  positions,  at  the  lowest 
temperatures,  lie  approximately  10-20%  below  the  predicted 
curve.  A  possible  explanation  of  the  observed  discre¬ 
pancies  will  be  given  below. 

A  magnetic  field  oriented  perpendicularly  to  the 

direction  of  the  tunneling  current  reduced  the  conduc¬ 
tance  anomaly.  This  effect  was  explained  on  the  basis 
of  the  Lorentz  force  induced  drift  of  the  electrons 
with  which  we  associated  a  relaxation  time  proportional 
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to  the  strength  of  the  applied  field.  No  observable 
changes  in  the  anomaly  were  observed  if  the  magnetic 
field  was  applied  parallel  to  the  direction  of  the 
current  flow. 

The  behaviour  of  the  conductance  with  the 
magnetic  field  is  instrumental  in  eliminating  the 
structure  around  zero  bias  as  being  due  to  paramagnetic 
impurities.  It  is  well  known  that  these  can  produce 
either  large  conductance  reductions  of  the  order  of 
100%  and  100  meV  wide  (Rowell,  Shen,  1966)  or  conduc¬ 
tance  peaks  (Wyatt,  1964)  depending  on  the  amount  of 
the  deposit  and  its  position  with  respect  to  the  barrier 
(Lythall ,  Wyatt,  1968;  El-Semary,  Rogers,  1972;  Appelbaum, 
Brinkman,  1970). 

Junctions  with  low  impurity  concentrations  produce 
conductance  peaks  which  are  temperature  and  magnetic 
field  dependent  (Shen,  Rowell,  1968;  Lythall,  Wyatt, 

1968) .  The  size  of  the  anomaly  is  usually  of  the  order 
of  10%,  being  a  few  meV  wide.  With  the  application  of 
the  magnetic  field,  the  conductance  peak  first  decreases 
and  with  higher  fields  splits  into  two  peaks,  the 
separation  of  which  is  proportional  to  the  field.  These 
results  have  been  explained  (Anderson,  1966;  Appelbaum, 
1966,  1967)  by  considering  the  interaction  between  the 
tunneling  electrons  and  the  magnetic  impurities.  Increas¬ 
ing  the  amount  of  the  deposit  leads  to  a  large  conductance 
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reduction  which,  however,  is  magnetic  field  independent 
(Rowell,  Shen,  1968).  It  is  now  believed  (Wolf, 
1974)  that  this  effect  is  more  likely  connected  with 
the  Gi aver- Zeller  (Zeller,  Giaver,  1969)  tunneling, 
rather  than  being  a  manifestation  of  the  interaction 
between  the  tunneling  electrons  and  the  interacting 
magnetic  impurities  (Gupta,  Upadhyaya,  1971) .  The 
dominant  current  flow  in  Giaver-Zeller  tunneling  is 
through  a  two  step  tunneling  process  in  which  an 
electron  first  tunnels  onto  the  particle  embedded 
inside  the  barrier  and  subsequently  tunnels  into  the 
second  electrode.  Since  the  localization  of  the 
electron  on  the  impurity  requires  an  activation  energy, 
this  process  leads  to  a  strongly  non  ohmic  behaviour. 

At  finite  temperatures  the  zero  bias  conductivity, 
a(0,T)  is  proportional  to  the  number  of  thermally 
excited  electrons  residing  on  the  particle.  This  is 
given  by 


a (0 ,T ) 


-b/kBT 
e  v 


0 


(4.10) 


where  E  is  the  activation  energy  of  the  excitation. 

At  this  point  we  call  attention  to  fig.  2.6, 
where  we  notice  that  the  zero  bias  blocking  conductance 
within  a  certain  temperature  limit  (depending  on  D  is 
a  linear  function  of  the  temperature.  Thus  the  behaviour 
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of  the  small  conductance  decreases  at  zero  bias  which 
is  linear  with  temperature  can  be  misleading  in 
believing  that  the  tunneling  is  through  particles 
(Josefowicz ,  Smith,  1973)  .  A  comparison  between  the 
observed  magnitudes  of  the  corresponding  anomalies 
(^100%  for  the  two  step  tunneling  vs.  <1%  for  the  non¬ 
equilibrium  effect)  and  their  magnetic  field  dependence 
helps  to  resolve  this  ambiguity. 

Let  us  turn  our  attention  to  the  possible  sources 
of  observed  differences  between  the  theory  and  the 
experiment.  One  of  these  arises  from  the  fact  that 
the  energy  and  temperature  dependent  electron-phonon 
relaxation  time  was  assumed  to  be  momentum  independent. 
There  are  several  sources  of  anisotropy.  The  most 
obvious  is  the  deviation  of  the  Fermi  surface  sphere. 
Another  results  from  the  anisotropy  in  the  electron- 
phonon  interaction  and  in  the  phonon  spectrum.  Recent 
calculations  (Carbotte,  Truant,  1974)  have  shown  that 
this  anisotropy  affects  strongly  the  value  of  the 
electron-phonon  relaxation  time.  Large  variations  over 
the  Fermi  surface  were  found  even  at  T  =  30  K.  For 

aluminum  the  ratio  of  Tep  ( 0 /■ $ *T ) /Tep  ( 0 ' 0 'T )  can  be  aS 
high  as  10  (K  E  (kp  ,  0  ,  <j> )  )  .  The  anisotropy  in  the  electron- 
phonon  relaxation  time  suggests  a  further  set  of  experi¬ 
ments  in  which  the  electrons  would  tunnel  into  single 


. 


124 


crystal  electrodes  having  a  definite  crystalographic 
orientation.  We  would  expect  a  certain  degree  of 
correlation  between  the  size  of  the  zero  bias  anomaly 
and  the  orientation  of  the  crystal  with  respect  to 
the  incoming  electrons. 

A  more  transparent  source  of  a  discrepancy 
arises  because  of  the  assumption  that  phonons  created 
by  the  relaxing  electrons  are  in  equilibrium  with  the 
thermally  excited  phonon  gas.  The  thermalization  of 
the  excess  phonons  was  assumed  to  be  immediate  with  the 
aid  of  a  perfect  heat  bath.  Thus  their  further  role 
in  the  deblocking  process  has  already  been  described 
by  the  previously  considered  electron-phonon  scatter¬ 
ing  process,  leading  to  a  relaxation  time  given  by 
eq.  (2.10).  Such  a  simplification  is  not  always  valid. 
If  the  lifetimes  of  the  emitted  phonons  were  "long 
enough"  an  incoming  electron  could  absorb  one  of  these 
and  thus  be  scattered  into  a  higher  energy  level.  This 
results  in  a  removal  of  the  electron  from  the  forward 
direction,  thereby  deblocking  the  occupied  state.  The 
process  is  in  sharp  contrast  with  the  deblocking  con¬ 
sidered  by  the  theory  previously  and  being  only  due  to 
the  interaction  between  the  tunneling  electrons  and 
thermal  phonons. 

It  is  illustrative  to  consider  the  order  of 
magnitudes  involved  in  order  to  obtain  an  e  stimate 
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whether  an  additional  deblocking  due  to  non-equilibrium 
phonon  absorption  is  possible.  The  rate  of  thermaliza- 
tion  through  the  surrounding  bath  can  be  characterized 
by 


L 

Tbath  v 

s 


(4.11) 


where  L  is  the  film  thickness  and  v  is  the  sound 

s 

velocity  in  the  metal.  For  a  1000  A  film,  and  v  = 

5  -11  . 

5xio  cm/p,  Tbath = 2  x 10  sec.  The  relaxation  time 

of  phonons  due  to  their  absorption  by  electrons  at 

T  =  0  K  is  given  by  equation  (4.12)  (Truell,  et  al,  1969) 


1 


T 


ph-e 


2_  Ef  ,  m  '  2 
9  it p  v  *2 

^O  S  T1 


(4.12) 


where  p  represents  the  density  of  the  metal.  At  1  meV 

1  5 

and  p  =  lOp/cm  ,  E.  =  10  eV,  v  =  5  x  10  cm 

O  J-  ^ 

t  ,  =  2  x 10  10  sec  . 

ph-e 


This  is  comparable  to  r^ath  an<^  the  relaxation  time  of 
electrons  due  to  their  interaction  with  the  equilibrium 
phonon  gas.  It  follows  that  electrons  can  be 

additionally  deblocked  by  scattering  with  the  non¬ 
equilibrium  phonon  gas.  Thus  in  steady  state  we  should 
solve  a  coupled  system  of  Boltzmann  transport  equations 
relating  the  non-equilibrium  situation  in  the  phonon 


and  electron  gas. 
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n  -  °  t  *  „ 

scatt  ph-e 


(4.13) 


0  =  (|f>  -  0 


I  3  f  \  I  /  3  f  >  ,  /  3  f  \ 

cMt  n.  +  (3t}  .  +  (3t} 

tunneling  e-ph  scatt 


(4.14) 


The  interpretation  of  different  terms  in  eqs.  (4.13)  and 
(4.14)  is  obvious. 

The  solution  of  the  coupled  equations  is  beyond 
the  scope  of  the  present  work.  Instead  of  this,  let 
us  see  what  possible  results  can  one  deduce  on  the  basis 
of  the  previous  considerations.  First  of  all  we  can 
crudely  estimate  the  importance  of  the  additional 
deblocking  mechanism.  Let  us  assume  that  the  effective 
relaxation  time  of  electrons  due  to  their  interaction 
with  phonons  can  be  written  as 


1 

T 


% 


TNo 


1 

T6N 


(4.15) 


where  N  =  No  +  6N .  -  t  -  gives  the  part  of  the 

lifetime  relaxation  which  is  associated  with  the  equi¬ 
librium  phonon  distribution  (eq.  2.10).  Setting 

T  =  T  (eq.  4.12)  we  further  assume  that  the  non- 

6N  ph-e  M 

equilibrium  phonons  relax  only  by  a  scattering  with  the 
transmitted  electrons.  The  amount  of  relaxation  due  to 
excess  of  phonons  is  negligible  only  if  >>  ^No* 

T  =  0  K,  recalling  eqs.  (2.12),  (4.12),  this  is  satisfied 
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(4.16) 


£  >>  .3  meV 


for  the  typically  used  quantities.  Practically  then  at 
v3  meV,  the  contribution  of  non-equilibrium  phonons  to 
deblocking  is  minimal. 

Secondly  at  T  =  0  K  we  expect  that  the  size  of 
the  conductance  dip  at  zero  bias  is  unaffected  by  the 
incorporation  of  the  non-equilibrium  phonons  into  the 
theory.  This  follows  from  the  fact  that  at  zero  bias 
there  are  no  non-equilibrium  phonons  created  since 
the  tunneling  current  is  zero.  As  voltage  is  imposed 
across  the  junction  the  zero  bias  anomaly  decreases 
sharper  because  of  the  additional  deblocking  mechanism 
previously  neglected  by  the  theory.  This  leads  to  a 
shift  of  V  ,  the  peak  position  in  da  /dV,  to  a  lower 
bias  with  a  simultaneous  increase  in  the  peak  amplitude 
With  increasing  temperatures,  the  influence  of  the  non¬ 
equilibrium  phonons  is  minimized  by  a  greater  thermal 
phonon  population. 

All  these  results  point  to  correction  factors 
which  when  included  in  the  theory  could  explain  the 
observed  experimental  discrepancies.  Since  the  nature 
of  this  discussion  is  purely  qualitative,  it  is  impos¬ 
sible  to  present  exact  values  of  these. 


CHAPTER  5 


RESONANT  ELECTRON  TUNNELING 

5 . 1  Introduction 

The  measurements  of  zero  bias  anomalies  on  some 
junctions  of  Mg/oxide /Mg,  Mg/oxide /Au  and  Mg/nit ride/Au 
have  produced  characteristics  which  had  marked  diffe¬ 
rences  from  those  shown  in  the  previous  chapter. 

As  a  representative  of  junctions  of  this  unusual 
behaviour  we  show  a  temperature  dependence  of  ™ 
for  junction  877A  (fig. 5.1).  Note  the  usual  set  of 
peaks  close  to  zero  bias  reflecting  the  discussed 
blocking  effect  with  very  weak  additional  structures 
seen  at  -3.5  meV  and  +1.8  meV.  The  structure  as  com¬ 
pared  to  the  background  is  less  than  5%,  with  the  back- 

-3  -1 

ground  itself  being  of  the  order  -5x10  meV  .  The 
order  of  magnitudes  involved  in  observation  of  such 
effects  again  documents  the  necessity  of  measurements 
at  lowest  attainable  temperatures  and  modulations  to  avoid 
thermal  and  instrumentational  smearing.  The  calculated 
relaxation  times  for  the  junction  877A  agree  well  with  the 
relaxation  times  for  samples  without  the  additional  struc¬ 
tures  . 

The  set  of  unexpected  peaks  in  prompted  an 

investigation  as  to  their  physical  origin  and  resulted 
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Fig .  5.1 


Illustration  of  a  junction  having  weak 
additional  structure  in  the  derivative 


of  the  conductance. 
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in  the  present  chapter.  To  substantiate  our  possible 
explanation  of  the  observed  effect  we  find  it  appro¬ 
priate  to  present  first  the  results  of  the  experiments 
on  junctions  showing  the  new  oscillatory  structure. 


5 . 2  Experimental  investigations  of  oscillations 

The  occurrence  of  the  new  structure  was  not 

restricted  to  junctions  of  Mg/Au  type.  In  fact,  we 

have  observed  the  additional  peaks  in  —  in  Mg/oxide/Mg 

and  Mg/oxide/Pb  junctions.  In  fig.  5.2  we  present  a 
du 

characteristics  for  Mg/oxide/Mg  junction  with  peaks 
which  are  better  resolved  and  are  more  numerous  in  the 
same  voltage  range  as  compared  to  fig.  5.1.  The 
corresponding  conductance  curve  is  shown  on  the  next 
fig.  5.3  with  noticeable  changes  in  its  shape  close  to 


da 

dV 


peak  positions.  It  is  interesting  to  note  that  in  one 


of  our  earlier  measurements  on  Mg/oxide/Pb  junctions  we 
have  seen  a  structure  similar  to  the  presently  investiga¬ 
ted  one  which  at  the  time  of  measurement  were  considered 
to  be  due  to  the  noise  in  the  junction. 

The  appearance  of  a  structure  in  Mg/Au  and  Mg/Pb 
junctions  gave  rise  to  a  suspicion  that  it  is  connected 
either  with  the  Mg  electrode  or  with  its  insulating  oxide 
In  view  of  this  fact,  we  have  decided  to  investigate  the 
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Fig.  5.2 


Conductance  vs  energy  for  a  Mg/Mg 
junction  having  oscillations. 
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Fig .  5.3 


Derivative  of  the  conductance  for 
a  Mg/Mg  having  oscillations 


883  A  Mg/Mg 
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new  effect  on  Mg/Au  junctions  since  the  likelihood  of 
their  successful  preparation  proved  to  be  the  greatest. 

We  have  accepted  the  observation  of  phonon  emission 
either  in  Mg  electrode  (Adler,  1969)  or  in  MgO  (when 
applicable)  (Adler,  1969)  as  a  sufficient  proof  that 
the  measured  current  was  due  to  tunneling.  Since 
the  (oeaks  in  da/dV  tend  to  have  an  oscillatory 
character,  we  shall  henceforth  refer  to  them  as 
"oscillations" . 

The  existence  of  oscillations  in  tunnel  junctions 
was  independent  of  the  techniques  of  sample  preparation 
as  described  in  Chapter  3,  and  the  occurrence  of  the 
additional  structure  even  in  junctions  prepared  under 
identical  experimental  conditions  was  unpredictable. 

We  have  considered  the  thickness  of  the  electrodes, 
junction  area,  the  chemical  composition  of  the  barriers, 
evaporation  rates,  glow  discharge  currents,  oxidation 
(nitrification)  times  and  pressures  as  variables  which 
could  have  influenced  the  final  properties  of  the  tunnel 
junctions.  The  magnesium  film  thickness  ranged  between 

O 

700-30000  A  with  the  gold  counterelectrode  of  600-3000  A 

2 

and  junction  areas  of  .01-. 05  mm  .  The  evaporation 
rates  varied  from  10  A/sec  to  1000  A/sec  at  pressures 
of  2xl0-6  atm.  The  combined  effects  of  the  glow  discharge 
current,  oxidation  (nitrification)  time  and  pressure  are 
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reflected  on  the  tunnel  junction  resistance  which  in 

4 

our  samples  ranged  between  15-8000  5!  at  He  tempera¬ 
tures.  The  decision  to  make  junctions  with  nitride 
barriers  was  taken  with  a  hope  of  altering  the  barrier 
composition  and  thus  possibly  changing  the  physical 
properties  of  the  metal-barrier  interface  which  we 
thought  to  be  responsible  for  the  appearance  of  the 
oscillations . 

In  fig.  5.4  to  fig.  5.12  we  present  some  typical 

characteristics  we  have  obtained  while  measuring  the 

properties  of  Mg/Au  junctions  having  magnesium  oxide 

or  nitride  as  their  insulating  barriers. 

To  illustrate  the  effect  of  the  temperature 

smearing,  we  show  ~  —  measured  at  1  K  and  4  K  for 

°o  av 

a  junction  with  weak  oscillations  (fig.  5.4).  First 

we  notice  a  decrease  in  the  blocking  effect  due  to  the 

elevated  temperature  with  a  continuous  smearing  of 

the  resolved  peaks  (at  1  K)  into  the  smooth  and  wide 

background  curve  at  4  K.  The  slight  overshooting  of  the 

1  K  curve  at  -5  meV  with  respect  to  the  background  is 

due  to  the  oscillatory  peak  at  -3.6  meV .  The  total 

variation  of  the  conductivity  between  0-5  meV  at  1  K, 

5 

represents  only  4  parts  in  10  . 

The  most  drastic  oscillatory  behaviour  was 
observed  in  875B  Mg/oxide/Au.  The  4  K  curve  (fig.  5.5) 
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Fig .  5.4 


Illustration  of  the  effect  of  the  temperature 

. _  _  .  .  .  1  da 

on  the  size  of  the  oscillations  in  — 
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Fig .  5.5 

Illustration  of  a  junction  having  well 
resolved  oscillatory  peaks  in  the 
conductance.  The  effect  of  the  temperature 
smearing  is  shown  as  well. 
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Fig .  5.6 

The  conductance  resp.  derivative  of  the 
conductance  vs  energy  for  875B  Mg/Au  measured 
to  ±25  meV.  The  arrows  indicate  the  expected 
phonon  emission  threshold. 
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again  passes  smoothly  through  the  oscillations  which  are 
strongly  pronounced  at  the  lowest  temperature.  In 
contrast  to  the  previously  seen  kinks  in  fig.  5.3  the 
conductivity  curve  now  displays  a  set  of  humps  almost 
equally  spaced.  The  asymmetric  offset  at  zero  bias  is 
probably  due  to  the  oscillatory  background.  The 
regularly  well  resolved  oscillations  persisted  to 
almost  ±10  meV  beyond  which  they  started  to  diminish 
in  size  (fig.  5.6).  The  arrows  indicate  the  expected 
phonon  emission  threshold  in  Mg.  For  the  positive 
applied  bias  the  conductivity  curve  has  a  series  of 
steplike  increases  which  are  almost  smeared  into  the 
structureless  background  beyond  10  meV.  Since  the  slope 
of  the  curve  is  smaller  on  the  negative  side  the  struc¬ 
ture  still  retained  its  humplike  behaviour  as  seen  pre¬ 
viously  in  fig.  5.3. 

Finally,  for  an  illustration  we  present  a  junc¬ 
tion  having  a  nitride  as  its  barrier  showing  oscillations 
in  ±25  meV  region.  The  4.2  K  curve  is  again  superimposed 
to  show  the  effect  of  the  temperature  smearing  (figs. 

5.7,  5.8)  . 

An  important  result  which  contributed  to  a  possi¬ 
ble  explanation  of  the  tunneling  curves  was  obtained  by 
an  adsorption  of  hydrocarbon  molecules  inside  the 
insulating  layer  during  the  barrier  formation.  Junctions 
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Fig .  5.7 


Illustration  of  the  oscillations  in  a  Mg/N2/Au 
junction  measured  to  ±25  meV .  The  arrows 
Indicate  the  expected  phonon  emission  threshold. 
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Fig .  5.8 


Derivative  of  the  conductance  vs  energy 
for  a  Mg/N2/Au  junction. 
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so  prepared  were  subjected  to  oil  vapours  during  oxida¬ 
tion  (nitrification)  by  passing  the  oxygen  (nitrogen) 
gas  through  a  flask  containing  light  machine  oil  heated 
close  to  its  boiling  point.  The  presence  of  adsorbed 
molecules  of  hydrocarbons  inside  the  barrier  drastically 
affected  the  tunneling  barrier  with  respect  to  the 
occurrence  of  the  oscillations.  In  all  junctions  of 
this  type  no  additional  structure  superimposed  on  the 
background  curve  was  seen. 

The  verification  of  the  oil  adsorption  was 
performed  by  an  observation  of  the  inelastic  spectrum 
of  the  junction  at  bias  range  where  previously  reported 
(Jacklevic  and  Lambe,  1968)  hydrocarbon  peaks  were 
located . 

To  gain  further  insight  into  the  nature  of  the 
oscillations  we  have  tried  to  alter  the  properties  of 
the  junctions  by  an  application  of  external  perturba¬ 
tions  such  as  annealing  the  junctions  to  room  temperatures, 
applying  voltages  up  to  1000  meV  across  the  junctions 
and  polarizing  the  interface  by  passing  DC  currents  along 
its  adjoining  electrode.  We  hoped  that  by  correlating 
the  newly  obtained  results  with  characteristics  measured 
before  the  application  of  the  external  influences  we 
could  extract  additional  information  which  would  help 
to  interprete  our  data.  Such  an  effort  proved  to  be 
quite  fruitful  as  it  will  be  seen  on  the  following  figures. 
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Junction  080B  when  measured  after  an  initial 
cooldown  to  1  K  had  oscillations  in  its  characteristics. 
These,  however,  disappeared  upon  annealing  the  junction 
to  room  temperature  while  still  keeping  it  in  a  vacuum 
(f igs . 5 . 9-5 . 10 )  .  The  annealing  produced  a  36%  decrease  in 
junction  resistance,  while  the  change  in  resistance 
over  the  spanned  voltage  region  approximately  remained 
the  same  (1.93%  and  2.02%).  The  sharper  conductivity 
increase  around  zero  bias  reflects  the  tail  of  the  first 
smeared  oscillation.  In  contrast  to  this  result,  an 
application  of  700  meV  (Mg  positive)  across  another 
junction  for  1  hour  produced  a  stronger  structure  which 
in  the  earlier  measurement  was  weakly  visible  (fig.  5.11). 
Besides  the  better  resolved  structure,  there  is  a  slight 
shift  in  peak  positions  (fig.  (5.12).  The  relative 
change  of  the  conductance  at  zero  bias  (o ^  2^)  “ 

^(0))/a^  1  (0)  =  0.0145  due  to  zero  bias  anomaly  in 
both  cases  remained  the  same  implying  that  the  boundary 
and  impurity  scattering  mechanisms  which  helped  to 
deblock  the  tunneling  electrons  in  the  metal  electrodes 
were  not  altered.  Thus  only  a  change  in  electron  transi¬ 
tion  rate  l/x^  is  responsible  for  the  change  in  the 
junction  resistance  (from  210.97fi  to  201.86  at  1  K) . 

Similar  results  were  obtained  by  passing  a  DC 
current  along  the  bottom  (Mg)  or  top  (Au)  electrode. 
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Fig .  5.9 

Comparison  of  the  conductance  curves  for 
0 8 OB  Mg/Au  before  and  after  annealing. 
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Fig .  5.10 


Comparison  of  the  derivative  of  conductance 
curves  for  080B  Mg/Au  before  and  after 


annealing . 
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Fig .  5.11 


Comparison  of  the  conductance  curves  for 
J107B  Mg/N2/Au  before  and  after  applying 
high  bias  across  the  junction. 
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Fig.  5.12 

Comparison  of  the  derivative  of  the 
conductance  curves  for  J107B  Mg/N2/Au 
before  and  after  applying  high  bias 
across  the  junction. 
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With  increasing  currents  the  resistance  of  the  junctions 
continuously  increased.  We  interprete  this  as  being 
due  to  the  burning  of  thin  edges  across  the  junction 
area  (Adler,  Kreuzer,  1972).  Simultaneously  with  this 
effect  we  have  seen  a  change  in  the  oscillation  pattern 
of  the  junctions. 

The  information  revealed  from  the  above  quoted 
results  represents  the  main  line  along  which  the  analy¬ 
sis  and  the  interpretation  of  the  data  will  be 
performed.  In  the  next  section  we  shall  consider  this 
and  suggest  a  physical  model  which  can  explain  the 
appearance  of  the  new  structure. 

5 . 3  Analysis  of  the  experimental  data 

Through  the  results  of  the  experimental  investi¬ 
gations  on  Mg/Au  and  Mg/Mg  junctions  a  new  electron 
transfer  mechanism  was  established  which  enhanced  the 
conductivity  of  the  tunnel  junctions.  The  existence 
of  the  associated  structure  was  found  to  be  independent 
of  the  thickness  of  the  electrodes  and  of  the  barrier. 
The  most  probable  bias  region  of  the  occurrence  of  the 
structure  was  found  to  be  within  25  meV  of  the  zero 
bias.  In  some  junctions,  however,  the  oscillations 
persisted  up  to  ±100  meV.  Since  the  amplitudes  of  the 
oscillations  rapidly  decreased  with  increasing  temper¬ 
atures  the  measurements  had  to  be  done  at  the  lowest 
attainable  temperatures. 
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The  appearance  of  oscillations  in  da/dV  with 
an  apparent  periodicity  of  0.8-2. 7  meV  led  to  an 
early  speculation  that  the  observed  effect  was  due 
to  resonances  associated  with  the  presence  of  the 
electron  standing  waves  in  Mg  (Jacklevic  et  al,  1972) . 
According  to  Jacklevic  et  al,  the  experimental  obser¬ 
vation  of  the  box  quantization  requires  a  preparation 
of  films  with  a  very  high  purity  and  smoothness.  This 
limited  the  thickness  of  Mg  films  in  their  experiments 
to  150  A-  500  A.  To  promote  smoothness,  the  films  had 
to  be  deposited  onto  substrates  cooled  to  -100°C-  150°C, 
with  a  subsequent  warmup  to  room  temperature  to  achieve 
recrystallization.  In  contrast  to  this  procedure  our 
Mg  films,  700  &  -  30000  A  thick >  were  evaporated  onto 
substrates  held  at  room  temperature.  We  would  then 
expect  that  the  stringent  conditions  imposed  on  the 
observation  of  the  standing  waves  were  not  reached. 
Furthermore,  correlation  between  the  thickness  of  the 
film  and  the  peak  separations  in  da/dV  has  not  been 
successful.  For  instance,  film  thicknesses  of  25000  it - 
30000  A  are  needed  in  order  to  resolve  oscillations  with 
a  periodicity  of  ^1  meV.  This  thickness  range  is  well 
above  the  thickness  of  the  Mg  film  in  junction  875 
(700  A) ,  in  which  the  oscillations  were  the  most  pro¬ 
nounced  from  all  our  samples.  These  considerations 
then  point  to  a  conclusion  that  the  structure  in  our 
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junctions  could  not  be  connected  with  Jacklevic  type 
resonances  . 

In  order  to  gain  further  insight  into  the 
nature  of  oscillations,  the  tunnel  junctions  were 
subjected  to  the  influence  of  external  perturbations. 

As  indicated  in  the  last  section,  a  thermal  cycling, 
an  application  of  high  biases  across  the  junction  or 
a  flow  of  a  DC  current  along  the  electrodes  produced 
irreversible  changes  in  the  junctions  which  resulted  in 
the  new  tunneling  characteristics.  These  could  be 
repeatedly  traced  over,  unless  a  new  perturbation  was 
reapplied.  We  have  to  mention, though, that  the  aging 
effects  in  thin  films  associated  with  the  thermal 
cycling  often  led  to  a  destruction  of  the  tunnel 
junction.  The  same  was  true  of  the  applied  DC  bias, 
resulting  in  the  dielectric  breakdown  of  the  insulator. 
Nevertheless,  we  feel  that  the  measurements  performed 
on  different  samples  revealed  enough  information  for 
their  mutual  intercomparison. 

We  attribute  the  changes  in  the  tunnel  junctions 
as  being  due  to  modifications  of  the  tunneling  layer  or 
of  the  interface.  We  further  assume  that  the  structure 
of  the  metal  electrodes  has  not  been  altered.  This 
is  supported  by  an  observation  that  the  relative  size 
of  the  conductance  anomaly  between  1  K  and  4  K  did  not 
change,  even  though  the  oscillatory  pattern  in  the 
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tunnel  characteristics  did  so.  The  new  tunneling  curves 
may  either  reflect  a  change  in  the  barrier  profile 
arising  from  the  additional  growth  of  the  insulating 
layer  (Pollack,  Morris,  1964)  with  a  simultaneous  change 
in  the  tunneling  path  or  result  from  the  spatial 
redistribution  of  impurities  inside  the  barrier.  Due 
to  the  exponential  dependence  of  the  transmission  pro¬ 
bability  on  the  thickness  and  height  of  the  barrier 
the  first  two  effects  merely  alter  the  magnitude  of 
the  tunneling  current  which  in  turn  determines  the 
structureless  background  (Simmons,  1963;  Stratton,  1962). 

A  considerable  variation  in  the  detailed  behaviour 
of  the  junctions  even  for  samples  prepared  under  iden¬ 
tical  condition  suggests  that  the  structure  is  connected 
with  the  impurity  content  of  the  insulating  layer. 
Clearly,  the  inelastic  impurity  assisted  tunneling 
( Jacklevic ,  Lambe,  1968)  associated  with  the  vibronic 
excitation  of  the  impurity  has  to  be  excluded.  This 
follows  since  it  is  not  expected  that  the  inelastic 
conductivity  enhancement  threshold  could  be  changed  by 
annealing  or  by  an  applied  DC  bias  (Geiger  et  al ,  1969). 
The  presence  of  an  impurity  inside  the  tunneling  layer, 
however,  often  leads  to  a  possibility  of  a  resonant- 
two  step  tunneling  process.  This  will  happen  if  the 
impurity  provides  an  attractive  potential  which  may 
have  a  bound  state  within  the  barrier  (Giaver-Zeller , 
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1968;  Parker,  Mead,  1969;  Gadzuk,  1970;  Duke,  1972). 
Then  under  favourable  conditions  a  new  tunneling 
channel  can  open  up,  in  which  the  tunneling  electron 
goes  from  the  metal  to  the  adsorbed  impurity  where  it 
is  partially  localized  and  then  tunnels  through  the 
reduced  barrier  into  the  second  electrode.  The 
degree  of  electron  localization  is  very  important  in 
distinguishing  between  two  possible  processes.  An 
impurity  state  with  a  deep  and  narrow  attractive 
potential  can  lead  to  a  capture  of  the  tunneling 
electron  with  a  characteristic  lifetime  longer  than 
the  barrier  tunneling  time. 

A  subsequent  decay  of  the  bound  state  with  the 
transmission  of  the  electron  into  the  second  electrode 
is  essentially  equivalent  to  Giaver-Zeller  tunneling. 

A  uniform  distribution  of  bound  state  energies  then 
leads  to  a  large  temperature  dependent  resistance  peak. 
In  contrast, the  effect  on  the  impurity  potential  on 
the  tunneling  electron  can  be  such  that  it  will  allow 
only  for  a  resonant  buildup  of  the  electron  wave  func¬ 
tion  in  the  potential  well  with  localization  times 
shorter  than  the  tunneling  times.  This  will  lead  to 
an  increase  in  the  electron  probability  amplitude 
inside  the  forbidden  region,  resulting  in  a  greatly 
enhanced  barrier  transmission  coefficient.  We  believe 
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that  this  type  of  the  electron  resonant  transmission 
is  responsible  for  the  structure  resolved  in  our 
tunnel  characteristics.  It  is  then  appropriate 
to  consider  a  model  which  explicitly  shows  how  such 
resonances  can  affect  the  conductivity  of  the  tunnel 
junctions . 

5 . 4  Resonant  tunneling  through  impurity  states 

For  simplification  let  us  consider  a  tunneling 
process  in  which  the  insulating  layer  contains  only 
one  impurity,  having  an  allowed  energy  level  of  a 
magnitude  eVQ,  inside  the  barrier.  For  further 
simplification  let  us  assume  that  the  metal  electrodes 
are  at  0  K  and  that  the  barrier  potential  can  be 
represented  by  its  average  value  4  (fig.  5.13a). 

As  the  bias  is  established  between  the  two  electrodes 
the  Fermi  levels  separate  and  the  impurity  level 
shifts  linearly  with  the  electric  field.  Past  the 
threshold  voltage 

V  =  I V  I  +  V  •  £ 

'o'  £ 

(where  r  is  the  position  of  the  impurity,  £,  thickness 
of  the  barrier)  the  impurity  lies  above  the  Fermi  level 
of  the  right  hand  electrode.  In  this  case  a  new  tunnel¬ 
ing  channel  opens  up  in  which  electrons  can  tunnel  from 


Fig.  5.13 


Schematic  potential  for  a  junction  having 
a  single  impurity  in  the  middle  of  the 
barrier : 

a)  no  applied  voltage 

b)  sufficient  applied  voltage  so  that  the 
resonance  channel  is  opened. 
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left  to  the  impurity  and  then  from  the  impurity  to  the 
right  (fig. 5. 13b). 

Let  us  consider  what  happens  to  the  wave  func¬ 
tion  of  an  electron  with  energy  equivalent  to  the 
trapping  level.  The  incoming  wave  after  penetration 
decays  exponentially  inside  the  barrier.  As  it 
approaches  the  trap  it  resonates  with  a  large  ampli¬ 
tude  inside  the  potential  well.  After  a  partial  locali¬ 
zation  it  decays  again  exponentially  towards  the  other 
side  of  the  barrier.  By  moving  the  potential  well  to 
the  middle  of  the  barrier  the  amplitude  of  the  incoming 
wave  can  be  made  equal  to  that  of  the  outgoing  wave 
(Bohm,  1964).  For  electron  energies  different  from 
the  resonance  level  the  wave  function  decays  exponen¬ 
tially  in  the  usual  way  while  passing  through  the  barrier. 

Gadzuk's  (Gadzuk,  1970)  treatment  of  the  described 
model  gives  the  following  expression  for  the  ratio  of 
conductivities  with  and  without  the  trapping  center: 

°-3.  =  qfi  .  i  exp  ( £  m  <tk)  ■  6  ( I V  j  -  2 1 V  I )  (5.1) 

°N  As  2  V  f. 

where  it  was  assumed  that  the  resonance  occurs  in  the 

low  voltage  region  (eV  <<  4>)  and  that  the  impurity  lies 

in  the  mid-plane  of  the  barrier.  INF  represents  the 

number  of  impurities  with  an  effective  cross  section  A^ . 

A  gives  the  junction  area.  0  (x)  is  a  usual  step 
s 
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function.  Setting  N.A./A  =  1 ,  4>  =  1  eV ,  £  =  30  A 

1  1  s 

the  ratio  of  a^/a^  at  the  threshold  of  impurity 

1 6 

assisted  tunneling  is  a^/a^  =  2  x io  .  However,  in 
real  situation  this  favourable  ratio  is  reduced 
considerably  by  smallness  of  N-^/A  . 

The  calculation  of  the  resonant  conductivity 
with  a  single  impurity  is  a  great  simplification  of 
the  tunneling  problem.  In  practice  the  impurities  are 
distributed  over  the  whole  volume  of  the  insulating 
layer.  It  is  then  reasonable  to  expect  that  the  re¬ 
sonance  levels  will  be  smeared  out  and  the  sharp 
conductivity  increase  will  be  further  reduced. 

Gadzuk ' s  calculation  shows  that  the  smearing  effect 
indeed  reduces  the  conductivity  increase.  Fig.  5.14a 
illustrates  the  ratio  of  the  distributed  impurity 
conductance  to  the  clean  junction  as  a  function  of 
applied  voltage,  treating  the  impurity  level  parame¬ 
trically  (N^A^/Ag  =  1).  Fig.  5.14b  shows  the  effect 
of  the  barrier  width  on  aR/a  .  The  broadening  of 
o^/ct,.  is  due  to  the  presence  of  a  finite  electric 
field  in  the  barrier  which  shifts  the  bound  state 
energy  of  an  impurity  by  an  amount  of  the  order  of  the 
voltage  drop  at  the  corresponding  impurity  site.  To 
avoid  complete  blurring  of  the  resonance  effect  it  is 
important  that  the  impurity  has  an  energy  level  close 
to  the  Fermi  level  and  that  it  is  located  in  a  few 
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Fig .  5.14 


The  ratio  of  distributed  impurity  conductance 
to  the  clean  junction  conductance  as  a  function 
of  voltage: 

a)  impurity  level  treated  parametrically 

b)  thickness  of  the  barrier  treated 
parametrically . 
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neighbouring  atomic  planes  inside  the  barrier.  Im¬ 
purity  clusters  lying  near  the  interface  yield  a 
contribution  which  is  very  asymmetric  with  respect 
to  the  applied  bias.  This  follows  because  the  shift 
of  the  energy  level  with  the  applied  voltage  is 
polarity  dependent. 

We  should  note  that  the  presented  treatment  is 
one  dimensional  and  assumes  that  the  impurity  has  only 
one  resonance  level .  In  reality  the  attractive 
potential  associated  with  the  impurity  can  have  more 
than  one  allowed  level  with  a  result  that  the  tunnel¬ 
ing  characteristics  will  show  a  series  of  peaks  which 
may  or  may  not  interact  with  each  other  to  bring  in 
additional  smearing  effect.  A  three  dimensional 
treatment  by  Hurault  (Hurault,  1971)  of  the  resonant 
electron  tunneling  predicts  that  the  transmission 
coefficient  varies  as  a  Lorentzian  with  the  peak  at 
the  resonance  energy  level.  The  width  of  the  Lorentzian 
is  equal  to  the  width  of  the  virtual  bound  state 

I T  I  ^  - y  X  |T.  I2 

(e  -  ed) 2  +  r2  imp 

R 


where  ER  is  the  resonant  level,  r  -  width  of  the  level 


and 


T  . 
imp 


is  the  resonance  enhancement  factor  depen¬ 
dent  on  the  position  of  the  impurity  and  its  concen¬ 
tration.  A  calculation  by  Combescot  (Combescot,  1971) 
confirms  Hurault' s  result. 
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5 . 5  Discussion 

On  the  basis  of  the  experimental  evidence  pre¬ 
sented  in  section  5 . 3  we  have  concluded  that  the 
interaction  between  the  tunneling  electron  and  the 
trapped  impurity  is  responsible  for  the  appearance 
of  the  additional  structure  in  our  tunnel  junctions. 
Schmidlin  (Schmidlin,  1965)  analyzed  the  effect  of  the 
trapping  centers  on  the  barrier  profile  and  found 
that  throughout  a  low  voltage  region  a  positive  ion 
substantially  reduces  the  barrier  high  in  its  locality 
At  large  biases  the  effect  of  the  impurity  is 
"smothered"  out  if  the  ion  is  close  to  the  positive 
electrode.  The  local  decrease  of  the  barrier  height 
around  the  impurity  can  then  be  instrumental  in 
establishing  a  resonant  transmission  of  electrons 
across  the  barrier  via  intermediate  energy  levels 
existing  inside  the  barrier. 

With  these  considerations,  the  changes  produced 
in  the  general  behaviour  of  the  junctions  by  annealing 
or  by  applying  a  DC  bias  can  readily  be  explained  by 
allowing  the  charged  impurities  to  migrate  through  the 
insulator.  The  change  in  the  oscillatory  pattern 
while  a  current  was  passed  along  the  electrode  is 
attributed  to  a  redistribution  of  charges  existing 
at  the  metal-insulator  interface. 


>1  - 
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A  convincing  confirmation  of  the  resonant 
impurity  assisted  tunneling  arises  from  the  obser¬ 
vation  of  the  behaviour  of  the  localized  energy  levels 
under  the  influence  of  a  magnetic  field.  In  fig.  5.15 
we  show  the  second  derivative  curve  of  a  previously 
shown  junction  (fig.  5.9)  measured  at  0  kG  and  60  kG, 
with  a  field  oriented  parallel  to  the  tunneling  current. 

The  characteristics  were  obtained  after  the  junction 
was  annealed  at  room  temperature.  Even  though  the 
oscillations  when  compared  to  the  first  measurement 
did  not  seem  to  exist,  a  10  fold  increase  in  sensitivity 
resolved  the  weak  structure.  The  prominent  feature  of 
the  60  kG  curve  in  fig.  5.15  is  the  appearance  of  two  peaks 
at  ^5  meV  as  a  result  of  the  spin  splitting  of  the 
previously  degenerate  atomic  level  (0  kG) .  The  effec¬ 
tive  g  factor  obtained  from  the  separation  is  g  =  2.75. 
With  the  exception  of  the  peak  at  ^16  meV  which  is 
identified  with  inelastic  phonon  emission  (Klein  et  al, 
1973)  the  other  peaks  appear  to  be  shifted  or  broadened 
under  the  influence  of  the  magnetic  field.  A  small 
energy  difference  in  the  separation  of  the  peak  at 
5  meV  suggests  that  the  strength  of  the  magnetic  field 
was  not  sufficient  to  achieve  complete  spitting  of  the 
other  atomic  levels.  Esaki  et  al  (Esaki  et  al,  1968) 
reported  a  similar  observation  of  localized  energy 
levels  in  the  forbidden  gap  of  PbTe .  Oscillations 
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Fig .  5.15 


Illustration  of  the  spin  splitting  of  an 
impurity  level. 


175 


ENERGY  (meV) 


176 


in  the  tunnel  characteristics  of  Pb/O/Pb  were  reported 
as  well  (Yanson  et  al,  1971).  These  were  directly 
correlated  with  the  impurity  contamination  of  the 
barrier  by  a  layer  of  uracil.  It  was  assumed  that 
tunneling  through  localized  energy  levels  was  res¬ 
ponsible  for  the  occurrence  of  the  oscillations. 

The  discussion  in  the  preceding  section  shows 

that  the  effect  of  the  impurities  is  most  pronounced 

if  they  are  distributed  in  a  few  atomic  planes  which 

are  well  separated  within  the  depth  of  the  insulator. 

The  sharpness  and  the  separation  of  the  experimentally 

observed  peaks  then  leads  to  a  conclusion  that  the 

impurities  do  not  interact  strongly  with  each  other. 

The  lifetimes  broadening  of  the  atomic  levels  of  the 

order  of  1  meV  gives  an  estimate  on  the  localization 

time  of  an  electron  inside  the  potential  well.  The 

-13 

time  on  the  order  of  6  x 10  sec  is  comparable  to 
the  accepted  magnitudes  of  the  barrier  transmission 
times  (Gadzuk,  Lucas,  1973;  Thornber  et  al,  1967). 
Recall  that  we  have "ab  initio"  required  such  a 
relationship.  The  bias  region  of  the  occurrence  of 
the  oscillations  confirms  the  theoretical  prediction 
that  their  observation  is  more  likely  if  the  inter¬ 
mediate  energy  levels  are  situated  in  the  vicinity  of 


the  Fermi  level. 


. 
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Even  though  the  origin  of  the  impurities  has 
not  been  established  a  few  considerations  may  help 
to  elucidate  this  question.  The  observation  of  a 
structure  in  Mg  junctions  irregardless  of  the  top 
electrode  and  of  the  barrier  leads  to  two  possibilities 
or  a  combination  of  both.  Into  the  first  category 
we  classify  a  contamination  of  the  barrier  (while  it 
was  grown)  by  foreign  atoms.  This  includes  the  un¬ 
avoidable  adsorption  of  residual  gases  inside  the 
vacuum  system.  Shen  (Shen,  1968)  suggested  that  various 
impurities  (A1 ,  Ni,  Pt)  can  be  introduced  into  the 
insulating  layer  if  the  glow  discharge  method  is  used 
for  the  barrier  formation.  The  contamination  was  caused 
by  materials  sputtered  from  the  electrodes.  The  Sn- 
Pb  and  Sn-Sn  junctions  showed  several  sharp,  peaklike 
1%  structures  in  the  conductance  of  the  tunnel  junctions 
which  he  identified  with  randomly  distributed  energy 
states  in  the  barrier. 

The  second  possibility  leads  to  an  assumption 
that  the  insulating  barrier  intrinsically  contains 
trapping  centers  arising  from  its  nonstoichiometric 
growth  (Fisher-Giaver ,  1961;  Pollack,  Morris,  1964; 
Esaki,  et  al ,  1968;  Sasaki,  1960).  Stevenson  and 
Hensley  (Stevenson,  Hensley,  1961)  concluded  that 
MgO  tends  to  exist  in  the  form  of  an  n-type  semi¬ 
conductor.  The  donors  responsible  for  this  behaviour 


. 


were  similarly  associated  with  the  nonstoichiometric 
composition  of  the  oxide.  We  suggest  further 
experiments  in  which  the  procedure  of  sample  preparation 
would  be  altered  in  order  to  eliminate  the  suspected 
reasons  behind  the  contamination  of  junctions.  We, 
similarly,  leave  open  the  answer  for  the  disappearance 
of  the  oscillations  with  the  oil  doping  of  the  junc¬ 
tions.  We  suggest  either  a  binding  of  uncompensated 
Mg  ions  into  complex  organic  molecules  or  a  much  simple 
reason  ,  arising  from  preventing  the  sputtering  of 
materials  from  the  discharge  electrodes  by  formation 
of  a  thin  layer  of  oil  on  the  surface  of  the  electrodes 


CHAPTER  6 


TUNNELING  INTO  BISMUTH 

6 . 1  Quantum  size  effect  in  Bi 

There  has  been  a  considerable  interest  in  Bi, 
both  theoretical  and  experimental  (Fal ' kovski j , 196 8 ) . 

Its  interesting  properties  result  largely  from  its 
carriers  which  are  very  light  and  few  in  number. 

The  fermi  level  of  the  electrons  is  ^25  meV 
(holes  v  11  meV)  and  the  effective  mass  along  some 
crystal  direction  is  two  or  three  orders  of  magnitude 
smaller  than  the  free  electron  mass  (Fal ' kovski j , 196 8 ) . 
Consequently  the  electronic  de  Broglie  wavelength  can 

O  O 

be  as  large  as  300  A  or  400  A.  In  comparison  the 
de  Broglie  wavelength  in  metals  is  on  the  order  of 

O 

the  lattice  period  which  is  at  the  most  a  few  A. 

If  the  de  Broglie  wavelength  of  the  conduction 
electrons  becomes  comparable  to  the  physical  dimensions 
of  the  sample,  the  electronic  energy  spectrum  should 
become  discrete  (quantum  size  effect)  with  the  separation 
between  the  energy  levels  being  determined  by  the 
sample  thickness  (Tavger,  Demikhovski j ,  1969).  The 
main  features  of  size  quantization  can  be  derived  from 
a  model  in  which  the  electrons  and  the  holes  are 
regarded  as  independent  particles  moving  in  an  infinite 
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square  well  potential.  The  single  particle  wave 
function  for  a  film  of  thickness  L  can  then  be 


written  as 


z 


ipn 


k  k  = 
x  y 
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z  x  y 


V2 
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sin 


nrz 
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(6.1) 


where  0  <  z  <  L  ,  and  the  film  has  infinite  dimensions 

z 

in  the  y  and  z  directions.  The  corresponding  energy 
spectrum  is 
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Thus  while  in  a  bulk  crystal,  the  electron  states  are 
given  by  points  which  fill  the  Fermi  sea  quasi- 
continuously ;  in  a  thin  film  the  Fermi  surface  breaks 
up  into  allowed  subsheets  parallel  to  k  -k  plane. 

Each  subsheet  corresponds  to  a  discrete  value  of  the 
wave  vector  along  the  quantized  direction  of  the  film 
(Garcia,  et  al,  1972).  These  consequences  have  been 
investigated  with  regard  to  electrical  conductivity, 
optical  properties,  galvanomagnetic  phenomena,  para¬ 
magnetic  susceptibility,  superconductivity  and  tunneling 
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phenomeha  (Tavger,  Demikhovski j ,  1969;  Larson,  1971). 

One  of  the  possible  methods  of  direct  experi¬ 
mental  observation  of  the  quantization  of  the  transverse 
motion  energy  in  a  film  is  the  use  of  tunneling  of  the 
electrons  in  a  sandwich  made  of  the  investigated  films 
(Gogadze ,  Kulik,  1965;  Kummel,  1968) .  In  the  tunnel 
junction,  the  films  of  quasi-discrete  energy  spectrum 
are  separated  by  a  potential  barrier.  When  an  external 
field  is  turned  on,  the  energy  levels  of  one  film  shift 
relative  to  the  levels  of  the  other  film.  Assuming 
elastic  transmission , the  current  through  the  barrier 
will  arise  whenever  the  populated  levels  of  one  film 
lie  opposite  the  free  levels  of  the  other  film.  As  a 
result,  a  current-voltage  characteristic  in  the  form 
of  a  series  of  peaks  is  obtained.  The  period  of  oscil¬ 
lations  is  approximately 


(6.3) 


where  m*  is  the  effective  mass  of  the  carrier.  The 
upper  bound  on  the  amplitude  of  the  oscillations  is 
estimated  to  be 
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(Gogadze,  Kulik,  1965).  It  is  illustrative  to  estimate 
the  order  of  magnitudes  involved.  In  order  that  the 
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oscillations  would  be  observable  let  us  require  that 


Tosc 

j -  ^  10  . 

non 

This  depending  on  the  Fermi  energy  puts  an  upper  limit 
on  the  thickness  of  the  film.  Table  6.1  illustrates 
the  calculated  values  of  the  required  thickness  and 
correspondingly  obtained  periods  of  oscillations,  AE . 

Table  6 . 1 

Thickness  and  periodicity  of  oscillations 
Material  m*  E^ 


typical  metal  m 


Bi  electrons  10 


o 


-2 


1  eV 


Required 

thickness 


120  A 


m 


o 


Bi  holes 


0.7m 


25  meV  7400  A 


15  meV  1140  A 


Calculated 

AE 

9  8  meV 

2.50  meV 

1.50  meV 


m  is  the  mass  of  the  free  electrons, 
o 

The  thickness  range  in  Bi  is  considerably  relaxed 
because  of  the  small  Fermi  levels  and  light  masses  of 
the  carriers.  Thus  according  to  theoretical  predictions 
it  is  possible  to  observe  the  QSE  directly  by  tunneling. 
However,  it  is  obvious  that  this  favourable  prediction 
is  further  minimized  through  thermal  smearing  and  scatter¬ 
ing  of  the  carriers  by  impurities  and  diffuse  boundaries. 
It  is  then  necessary  to  keep  the  films  sufficiently 
homogeneous  and  pure.  Tunneling  into  single  crystals 
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or  films  containing  large  grains  of  well  defined 
orientation  therefore  increases  the  likelihood  of  the 
observation  of  the  QSE .  The  existence  of  quantized 
energy  levels  in  thin  metal  films  was  well  documented 
through  ingenious  interpretation  of  tunneling  results 
of  Jacklevic  et  al  (Jacklevic  et  al ,  1971).  They 
have  noticed  that  even  though  the  required  smoothness 
of  the  metal  films  is  not  experimentally  feasible, 
certain  energy  levels  (for  instance  at  the  band  edge) 
are  stationery  with  respect  to  thickness  perturbations. 
These  then  are  easily  determined  by  tunneling  experi¬ 
ments  . 

A  direct  observation  of  the  quasi-discrete  energy 
level  spectrum  by  tunneling  was  claimed  by  Russian 
experimentalists (Korneev  et  al,  1970;  Lutskii  et  al, 
1966).  Korneev  et  al  have  observed  oscillations  in 
Ag-Bi  structures  in  which  the  electrodes  were  condensed 
onto  two  separate  mica  substrates  and  pressed  together. 
The  measurements  were  performed  at  nitrogen  temperatures. 
The  thickness  of  Bi  films  was  greater  than  1000  A. 
Lutskii  et  al  reported  the  observation  of  the  QSE  in 
structures  in  which  the  thickness  of  the  Bi  films  ranged 
between  800-1300  A.  In  both  cases  the  steplike  increases 
in  the  tunnel  characteristics  allowed  a  successful  cor¬ 
relation  with  the  thickness  of  the  film. 

The  purpose  of  our  experiments  was  to  repro¬ 
duce  these  results  through  tunneling  experiments  in 
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systems  containing  Bi  films  in  which  the  insulating 
layer  was  formed  by  a  solid  dielectric. 

6 . 2  General  considerations 

The  investigations  were  performed  on  Al/Bi, 

Mg/Bi,  Bi/Bi  and  Bi/Pb  junctions.  They  were  pre¬ 
pared  according  to  the  procedure  described  earlier 
(section  3.2).  The  substrates  usually  were  heated 
to  80 °C  by  a  resistive  heater  in  order  to  promote 
formation  of  large  grains  in  the  evaporated  film. 

Even  though  Bi  films  evaporated  onto  glass  slides 
generally  are  polycrystalline,  the  trigonal  axes  of 
their  grains  are  usually  within  20%  of  the  film 
normal  (Komnik,  Bukhstab,  1968).  Condensation  of 
Bi  films  onto  heated  mica  creates  larger  crystallites, 
the  trigonal  axes  of  which  are  perpendicular  to  the 
substrate  (Komnik,  Pilipenko,  1971).  Hauser  and 
Testardi  (Hauser,  Testardi,  1968)  reported  similar 
observations.  Electron  microscope  study  showed  that 
the  lateral  sizes  of  the  microcrystals  were  on  the 
order  of  1500  A  -  3000  K.  The  extent  of  grain  forma¬ 
tion  was  not  very  critical  on  the  evaporation  rate. 

The  annealing  of  films  during  oxidation  did  not  increase 
the  grain  sizes  appreciably  (20%-30%)  . 
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The  thickness  of  Bi  electrodes  ranged  between 
1500  A-  15000  A.  Films  thinner  than  1500  A  after 
oxidation  had.  very  high  resistances  (>5  kft)  . 

Junctions  made  of  these  in  most  cases  were  not  stable 
and  were  easily  destroyed  by  biasing  them  only  to  a 
few  meV .  The  counterelectrodes  of  A1 ,  Mg,  Pb  were 
approximately  1000-  5000  A  thick. 

The  detection  of  an  energy  gap  in  Pb  by 
tunneling  gave  us  a  confidence  that  the  dominant 
mechanism  of  a  current  transport  across  the  oxide 
of  Bi  was  due  to  tunneling.  The  observation  of 
inelastic  phonon  assisted  tunneling  in  Bi/Bi 
samples  similarly  led  us  to  believe  that  by  oxidation 
of  Bi  in  oven  we  could  prepare  "good"  tunnel  junctions. 

6 . 3  Discussion  of  results 

Early  measurements  on  Bi/Bi  junctions  showed 
a  reproducible  oscillatory  type  structure  in  the 
vicinity  of  zero  bias.  We  illustrate  this  in 
figs.  6.1  and  6.2.  Fig.  6.1  shows  the  conductance 
curves  for  three  Bi/Bi  junctions.  The  Bi  films  in 
junctions  007  and  016  were  5000  A  thick.  The  films 
in  001  were  4000  A  thick.  The  corresponding  derivatives 
da/dV  are  in  the  next  figure  (fig.  6.2).  The  amplitudes 
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Fig .  6.1 


Conductance  vs  energy  for  three  Bi/Bi 
junctions.  The  thicknesses  of  films  were 
following 


001  -  4020  A/ 4020  A 
007  -  5360  A/5360  A 
016  -  5420  A/5360  A 
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da/dV  vs  energy  for  three  Bi/B 
junctions 
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of  oscillations  in  001,  in  comparison  to  the  two 
other  junctions,  are  small  and  barely  visible.  For 
positive  applied  voltages  (top  film  positive)  the 
oscillations  had  practically  the  same  periodicity. 

This,  however,  was  not  true  for  the  bias  voltages  of 
negative  polarity  (top  film  negative) .  The  explanation 
of  this  effect  is  probably  due  to  difference  in  the 
film  structure,  since  the  bottom  films  were  always 
subjected  to  an  annealing. 

Fig.  6.3  illustrates  dc/dV  for  007  in  a  greater 
detail.  The  general  features  of  oscillations  seen 
in  fig.  6.3  were  well  reproduced  even  in  junctions 
where  the  electrodes  were  15000  A  thick.  In  all 
samples  the  regular  oscillations  always  disappeared 
beyond  the  phonon  emission  threshold  at  3  meV  (Rowell 
et  al ,  1969).  Taking  into  account  this  reproducibility 
we  believe  that  the  structure  cannot  be  connected 
with  impurity  assisted  tunneling  discussed  in  Chapter  5. 
However,  since  the  periodicity  of  oscillations  did  not 
change  with  increased  thickness  of  Bi  films,  we  can 
immediately  conclude  that  the  structure  cannot  be 
connected  with  the  observation  of  quantum  size  effect. 
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Fig .  6.3 

da/dV  vs  energy  for  007  Bi/Bi  shown  separately 


ENERGY  (meV) 


193 


6 . 4  Phonons  in  Bi 

The  observation  of  phonon  emission  by  tunneling 
in  carefully  prepared  single  crystals  of  Bi  has  been 
reported  by  Esaki  et  al  (Esaki,  et  al,  1968a).  The 
conductance  characteristics  of  Bi/A^O^/Au  tunnel 
junctions,  in  which  Bi  was  oriented  along  its  trigonal 
axis,  showed  a  set  of  humps  located  approximately  at 
3±1  meV,  6±1  meV,  10±1  meV  and  14±1  meV,  regardless  of 
the  polarity  of  applied  voltage.  These  were  correlated 
with  previously  measured  phonon  energies  of  the  four 
propagation  modes  in  the  trigonal  direction  (Kotov 
et  al ) : 

transverse  acoustic  (TA)  =  4.8  meV  , 

longitudinal  acoustic  (LA)  =  7.4  meV  , 

transverse  optic  (TO)  =  12.6  meV  , 

and  longitudinal  optic  (L0)=  13.4  meV  , 

at  the  Brillouin  zone  boundary, 

and  TO  =  9.2  meV  and  LO  =  12.5  meV 

at  the  zone  center. 

Peaks  in  da  /dV  at  3  and  12  meV  reflecting  the  inelas- 
e 

tic  excitation  of  Bi  phonons  in  Al/Bi  junctions  were 
observed  by  Rowell  et  al  (Rowell,  et  al ,  1969)  . 

Results  similar  to  those  quoted  above  were 
obtained  in  our  samples  as  well.  However,  in  contrast 
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Fig .  6.4 

Observation  of  inelastic  phonon  assisted 
tunneling  in  a  Bi/Bi  junction  prepared  on 
mica.  The  numerals  show  the  location  of 
phonon  peaks. 
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Fig .  6.5 

Observation  of  inelastic  phonon  assisted 
tunneling  in  a  Mg/Bi  junction.  The 
numerals  show  the  location  of  phonon  peaks . 
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Fig.  6.6 

Observation  of  inelastic  phonon  assisted 
tunneling  in  an  Al/Bi  junction.  The 
numerals  show  the  location  of  phonon  peaks. 
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to  the  observation  of  Rowell  et  al,  we  were  able  to 
detect  phonons  which  energetically  lie  between  3  and 
12  meV. 

According  to  Esaki  et  al  junctions  with  poly¬ 
crystalline  Bi  films  did  not  exhibit  the  fine  phonon 
structure  in  their  tunnel  characteristics.  Further¬ 
more,  they  did  not  see  any  structure  in  junctions, 
the  resistances  of  which  were  below  5  kfi.  In  com¬ 
parison  the  resistances  of  our  samples  were  less  than 
5  kfi.  The  detection  of  phonons  in  our  polycrystalline 
Bi  films  therefore  complements  the  results  of  Esaki 
et  al . 

The  assignments  of  phonon  emission  thresholds 
according  to  Kotov  et  al  are  following  :— 


Identif 

Peak  number 
1 
2 

3 

4 


Table  6 
cation  of  phonon 

Position (meV) 
3.1 

4.8 

7.8 
12.3 


2 


peaks  in  Bi/Bi  junction 

Possible  Cryst. 

identification  direction 


LA 

TA 

LA 

TO  ,LO 


2  fold  axis 

3  fold  axis 
3  fold  axis 
3  fold  axis 


Similar  identification  is  performed  in  a  Mg/Bi  junction 
(fig .  6.5). 
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Table  6.3 


Identification  of  phonon  peaks  in  Mg/Bi 


Peak  number 

Position (meV) 

Possible 

identification 

cryst . 
direction 

1 

3.3 

LA  2 

fold 

axis 

2 

7.7 

LA  3 

fold 

axis 

3 

9.7 

TO  3 

fold 

axis 

4 

12.4 

TO , LO  3 

fold 

axis 

5 

18.4 

transverse  Mg 

6 

27.5 

longitudinal  Mg 

The  unmarked 

peak  is  due  to 

zero  bias  anomaly. 

Finally ,  we  present  the  results  of  tunneling  into  an 
Al/Bi  junction  (fig.  6.6). 

Table  6.4 

Identification  of  phonon  peaks  in  Al/Bi 

Possible  Cryst. 

Peak  number  Position (meV)  identification  direction 

LA 
TA 
LA? 

LO 


1 

2 

3 

4 


3.7 

6.5 

8.5 
11.4 


2  fold  axis 

3  fold  axis 
3  fold  axis 
2  fold  axis 


5 


15.1 


2  phonon  processes? 


- 


When  determining  the  exact  peak  positions  above 


a  small  correction  factor  (.18  meV)  due  to  super¬ 
conductivity  of  the  Al-electrode  was  subtracted 
from  the  experimentally  measured  peak  positions. 

The  peak  at  6.5  meV  is  identified  only  tentatively 
in  the  view  of  uncertainty  in  exact  determination  of 
maximum  in  do^/dV.  The  peak  at  8.5  meV  is  above  the  LA 
phonon  (7.3  meV)  emission  threshold.  Both  peaks,  how¬ 
ever,  at  4.2  K  merged  into  a  broad  shoulder  around 
8  meV  which  is  in  the  vicinity  of  the  excitation  energy 
of  the  LA  phonon.  The  common  feature  of  all  three 
samples  is  the  appearance  of  phonon  peaks  in  the 
vicinity  of  3  and  12  meV.  This  suggests  that  the 
coupling  between  electrons  and  phonons  is  strongest 
at  these  energies.  We  expect  that  the  resolution  of 
phonon  peaks  at  different  energies  depends  on  the  size 
and  orientation  of  grains  comprising  the  evaporated 
Bi  film. 


6 . 5  Other  interesting  effects 

In  connection  with  attempts  to  observe  QSE  in 
Bi  we  have  performed  tunneling  experiments  in  which 
the  junctions  were  biased  up  to  500  meV .  All  Al/Bi 
and  Mg/Bi  junctions  have  shown  the  characteristic  W- 
shape  (Fig.  6.7).  The  Bi/Bi  junctions  on  the  other 
hand  show  the  more  parabolic  shape  associated  with 
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Fig .  6.7 

Normalized  conductance  vs  energy  for  typical 
Al/Bi,  Mg/Bi  and  Bi/Bi  junctions.  The  Mg/Bi 
curve  has  been  displaced  upwards  by  0.2  for 
clarity.  Negative  bias  indicates  that  the 
bottom  layer  is  negative  with  respect  to  the 


counterelectrode . 
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symmetrical  tunnel  junctions.  We  feel  that  the  W- 
shape  arises  from  the  large  difference  in  Fermi 
energies  between  the  base  layer  and  Bi  counter¬ 
electrode  . 

In  the  intermediate  energy  range  (<100  meV) 
we  have  observed  a  structure  in  Bi/Bi,  Mg/Bi, 

Bi/Pb  and  Pb/Bi  junctions  which  in  certain  extent 
correlated  with  the  known  band  edges  of  Bi  deter¬ 
mined  by  other  observers  (Hauser,  Testardi,  1968; 
Esaki  et  al,  1965). 


CHAPTER  7 


SUMMARY  AND  CONCLUSIONS 

7 . 1  Discussions 

The  experimental  work  presented  in  this  thesis 
is  divided  into  three  chapters. 

In  Chapter  4  we  have  presented  the  results  of 
measurements  of  small  conductance  decreases  occurr¬ 
ing  around  zero  bias.  A  recent  theory  (Trof imenkof f 
et  al,  1972)  has  explained  these  through  non-equilibrium 
treatment  of  the  tunneling  process.  In  particular,  the 
related  structure  in  a  and  da/dV  have  been  accounted 
for  by  the  blocking  of  otherwise  available  electron 
states  due  to  finite  electron  relaxation  times  in  metal 
electrodes.  The  theory  assumes  that  the  deblocking  of 
occupied  states  is  through  scattering  of  electrons  via 
lattice  vibrations  and  impurities.  These  processes  are 
characterized  by  a  temperature  and  energy  dependent 
electron-phonon  relaxation  time  x  and  a  constant 
impurity  scattering  time,  x^  The  temperature  depen¬ 
dence  of  the  zero  bias  anomaly  is  described  only  by  one 
parameter  F,  the  knowledge  of  which  allows  one  to 
extract  the  impurity  scattering  time  x.^,  if  the  strength 
of  the  electron-phonon  interaction  is  known  beforehand. 
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We  have  performed  detailed  measurements  on  Mg/Mg, 

Mg/Au  and  Al/Al  junctions  and  found  that  the  theory 
agrees  well  with  experiment.  A  magnetic  field 
oriented  perpendicularly  to  the  direction  of  tunnel¬ 
ing  current  reduced  the  conductance  anomaly.  This 
effect  was  explained  on  the  basis  of  Lorentz  force 
induced  drift  of  tunneling  electrons  with  which  we 
have  associated  a  scattering  time,  Tdrift'  ProPor_ 
tional  to  the  strength  of  applied  field.  No 
observable  changes  in  the  anomaly  were  detected  if 
the  magnetic  field  was  applied  parallel  to 
the  tunneling  current. 

While  investigating  the  zero  bias  anomaly 
we  have  observed  additional  oscillatory  type  struc¬ 
ture  in  characteristics  of  some  tunnel  junctions. 

The  structure  was  found  to  exist  in  Mg/Mg,  Mg/Au  and 
Mg/Pb  type  junctions,  irregardless  of  electrode  thick¬ 
ness  and  barrier  preparation.  With  the  aid  of  exter¬ 
nal  perturbations  the  properties  of  junctions  could 
be  irreversibly  altered  resulting  in  new  characteristics. 
Adsorption  of  hydrocarbons  during  formation  of  barrier 
layer  prevented  the  appearance  of  oscillations.  From 
these  observations  we  have  concluded  that  a  contamina¬ 
tion  of  the  insulating  layer  by  charged  impurities  is 
responsible  for  the  observed  structure.  It  was  suggested 
that  the  impurities  could  alter  the  barrier  profile  in 
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such  a  way  as  to  allow  for  a  resonant  transmission  of 
tunneling  electrons.  The  effect  of  external  pertur¬ 
bations  is  then  understood  on  the  basis  of  rearrange¬ 
ment  of  impurities  inside  the  barrier.  Another 
possible  explanation  of  the  observed  effect  arises 
if  we  consider  a  creation  of  space  charge  trapping 
clusters,  which  could  be  instrumental  in  capture  of 
tunneling  electrons.  However,  until  the  exact  nature 
of  impurity  contamination  is  not  found,  the  results  of 
discussion  given  in  Chapter  5  will  remain  purely 
speculative . 

In  Chapter  6  we  have  briefly  commented  on  an 
unsuccessful  attempt  to  observe  quantized  energy 
levels  of  Bi  by  tunneling.  As  a  byproduct  of  mea¬ 
surements,  we  have  observed  the  excitation  of  phonons 
in  polycrystalline  Bi  films  with  Bi/Bi,  Mg/Bi  and 
Al/Bi  junctions.  These  results  augment  the  previously 
published  data  and  simultaneously  provide  better 
resolution  for  some  of  the  observed  peaks.  The  high 
energy  conductance  curves  (^500  meV)  of  Al/Bi  and 
Mg/Bi  junctions  displayed  the  usual  W-shape  associated 
with  such  junctions.  On  the  other  hand,  Bi/Bi  junctions 
showed  characteristics  which  are  typical  of  symmetrical 
junctions.  This  interesting  result  requires  further 
experiments  for  its  proper  interpretation. 
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7 . 2  Suggestions  for  further  work 

Concerning  the  zero  bias  anomaly  we  propose  to 
investigate  these  areas  as  possible  extensions  of  the 
presented  work. 

Evaluate  the  dependence  of  the  zero  bias  anomaly 
on  the  crystal  structure  of  the  metal  electrons.  In 
particular,  we  would  expect  a  correlation  between  the 
magnitude  of  the  conductance  dip  and  the  anisotropy 
in  the  electron-phonon  relaxation  time.  This,  of  course, 
requires  a  preparation  of  single  crystal  junctions.  We 
include  here  the  dependence  of  the  zero  bias  anomaly  on 
film  thickness.  The  experimental  results  gave  quali¬ 
tative  agreement  with  the  theory.  A  more  detailed 
investigation  should  possibly  extract  the  amount  of 
blocking  due  to  decreased  film  thickness.  In  connection 
with  this  we  suggest  backing  of  electrodes  with  a 
silver  layer  in  order  to  change  the  amount  of  specular 

scattering . 

A  further  study  should  show  the  effect  of  alloy 
ing  on  the  size  of  zero  bias  anomaly.  With  such  experi¬ 
ments  we  can  change  not  only  the  impurity  scattering 
time  xif  but  the  electron-phonon  coupling  strength  as 
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Finally,  as  it  has  already  been  shown 
(Wattamaniuk ,  1972)  the  non-equilibrium  aspect  of 
the  electron  tunneling  can  be  extended  to  include 
a  tunneling  into  superconducting  metals.  In  par¬ 
ticular,  for  the  case  of  two  identical  superconductors, 
the  continuous  increase  in  the  tunnel  current  just 
above  twice  the  gap  edge  (T  ^  0  K)  has  been  explained 
as  due  to  blocking  of  the  incoming  quasi-particles. 

It  should  be  noted  that  the  quasi-equilibrium  theory 
predicts  a  discontinuous  jump  which  is  experimentally 
not  observed.  Further  experiments  are  needed  to  clarify 
this  situation. 

Regarding  the  impurity  assisted  tunneling,  the 

objective  of  further  experiments  is  to  determine  the 

( 

origin  of  barrier  contamination.  This  requires  a 
different  method  of  barrier  formation;  thermally  grown 
barrier  eliminates  a  possibility  of  contamination 
through  sputtering  action  during  glow  discharge. 
Experiments  should  be  performed  in  order  to  determine 
a  mechanism  by  which  adsorption  of  hydrocarbons  in  the 
insulating  layer  prevents  the  occurrence  of  oscillations. 

Since  our  attempts  to  observe  the  quantized 
energy  levels  directly  by  tunneling  have  failed  we 
suggest  further  experiments  in  which  surface  roughness 
of  Bi  films  is  minimized  and  better  controlled.  This 
includes  tunneling  into  single  crystals  of  Bi. 
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Furthermore,  tunneling  in  Bi/oxide/metal  and  metal/ 
oxide/Bi  type  junctions  appears  to  be  a  unique  method 
of  proper  identification  of  a  structure  in  character¬ 
istics  of  these  junctions. 
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A  system  which  is  capable  of  resolving  changes  in  a  tunneling  junction  conductance  a  of  about  one  part  in  106  as 
well  as  providing  calibrated  values  of  dtr/dV  is  described.  Methods  of  simultaneous  analog  and  digital  recording  on 
magnetic  tape  suitable  for  subsequent  computer  analysis  are  given  .  The  flow  chart  for  the  computer  programming  is 
also  presented  along  with  a  typical  result. 


INTRODUCTION 

LECTRON  tunneling  experiments  have  been  very 
valuable  in  obtaining  detailed  information  about  the 
interactions  giving  rise  to  superconductivity.1'2  More  re¬ 
cently,  this  technique  has  been  applied  to  the  study  of 
inelastic  excitations  owing  to  phonons  and  molecular 
vibrations  in  metal-insulator-metal  (M-I-M)  tunnel 
junctions  in  which  both  metals  are  normal.2-7  The  quanti¬ 
ties  of  interest  in  these  studies  are  the  dynamic  conduct¬ 
ance  a  =  dI/dV  of  the  junction,  and  its  derivative  dcr/dV 
=  (PI/dV 2  as  a  function  of  junction  bias  voltage  V. 

This  paper  describes  instrumentation  capable  of  pro¬ 
viding  measurements  of  a  and  do/dV  and  recording  these 
on  magnetic  tape  suitable  for  subsequent  computer 
analysis.  The  system  is  capable  of  resolving  changes  in  cr 
of  about  one  part  in  105  as  well  as  providing  calibrated 
values  of  da/dV.  The  method  of  measurement  employs  a 
combination  of  harmonic  detection  and  bridge  techniques 
similar  to  that  described  earlier  by  Adler  and  Jackson8 
(henceforth  referred  to  as  I)  along  with  a  data  acquisition 
system  which  outputs  the  data  on  digital  magnetic  tape. 
The  method  of  calibrating  da/ dV  with  the  aid  of  a  digital 
computer  is  also  discussed. 

The  paper  is  divided  into  three  parts  and  an  Appendix. 
Following  this  introduction,  Sec  I  describes  the  circuit 
layout  and  discusses  the  bridge,  sweep,  and  modulation 
circuit  inasmuch  as  these  differ  from  those  described 
earlier  in  I.  The  digital  recording  formatting  and  buffering 
are  also  discussed.,  Section  II  deals  with  the  calibration  of 
a  and  da/dV  and  develops  an  integral  equation  scheme 
whereby  da/dV  may  be  calibrated  in  absolute  terms.  A  flow¬ 
chart  for  computer  programming  of  this  calibration  is  also 
given.  The  final  section  summarizes  the  operating  pro¬ 
cedure  and  gives  typical  results.  At  the  end  of  the  paper 
there  is  an  Appendix  which  deals  briefly  with  the  errors 
which  can  occur  owing  to  excessive  harmonic  distortion 
in  the  modulation  supply. 

I.  THE  CIRCUITRY 

The  analog  part  of  the  circuitry  described  in  this  section 
is  an  outgrowth  of  the  bridge  described  in  I  and  we  will 
consider  here  only  those  aspects  of  the  bridge  not  pre¬ 
viously  discussed.  Figure  1  shows  a  schematic  of  the  bridge 


enclosed  in  a  rectangle  of  dashed  lines.  This  bridge  differs 
from  the  earlier  version  in  the  following  ways. 

(1)  A  variable  capacitance  Cd  (100  nF  in  steps  of  0.1  nF) 
has  been  added  to  compensate  for  the  junction  and  stray 
capacitance. 

(2)  The  output  of  the  bridge  is  now  transformer  coupled 
to  the  preamplifier  (Triad  G-10  transformer).  This  trans¬ 
former  is  preceded  by  20  mI  polystyrene  capacitors  to 
prevent  the  possibility  of  direct  current  altering  the  trans¬ 
former  core  hysteretically. 

The  bridge  arms  R  are  chosen  so  that  R»Rr,  where  Rt 
is  the  junction  resistance.  The  sweep  and  modulation 
signal  circuits  are  showm  in  Fig.  2.  If  the  current  modula¬ 
tion  across  the  junction  is  given  by  b  cose at,  where  8  is  a 
small  constant  current  amplitude,  then  from  I  we  can 
write  the  voltage  output  across  the  transformer  at  angular 
frequency  to  as 

V„=(b/$)[Rd- (dV/dI)iJ  cosa d. . .,  (1) 

where  (dV/dI)it  is  the  dynamic  resistance  of  the  junction 
biased  by  a  current  /o •  The  factor  d  has  the  form  1  -\-Rd/R 
+2Rd/Z  where  Z  is  the  transformer  impedance.  The 
modulation  level  is  kept  low  so  that  the  instrumental 
smearing  is  of  the  order  of  the  thermal  smearing  (kT  at 
1  K  is  86  mV  rms).  Since  we  work  primarily  at  liquid  helium 
temperatures  the  signal  levels  are  very  low  so  the  bi  idge 
has  to  be  well  shielded  and  leads  kept  short.  The  second 
harmonic  at  the  transformer  will  be  given  by 

V ->a {b'2 (d2V /dP) i a  coslut.  ...  (2) 

Since  the  bridge  is  operated  keeping  a  constant  cunent 
through  the  junction  ( R»Rr )  the  measurements  are  re¬ 
stricted  to  junctions  where  V  is  single  valued  as  a  function 
of  current  I.  That  is,  negative  resistance  regions  such  as 
those  which  occur  in  M-I-M  junctions  in  which  both 
metals  are  superconducting  below  2  mV  cannot  be  handled 
by  this  circuity.  Such  negative  resistance  regions  can  how¬ 
ever  be  measured  by  other  methods.8'9  Errors  owing  to 
small  amounts  of  harmonic  distortion  in  the  modulation 
supply  can  be  serious  if  the  bridge  is  operated  far  off 
balance  (very  nonlinear  junctions) ;  these  are  discussed  in 
the  Appendix  and  the  above  bridge  equations  are  re¬ 
derived  considering  a  small  component  of  second  har¬ 
monic  present. 
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We  turn  our  attention  now  to  the  analog  and  digital 
recording  of  the  data.  The  bridge  signal  is  amplified  by  a 
PAR  type  113  preamplifier  followed  by  a  PAR  model  122 
lock-in  amplifier,  whose  dc  output  is  fed  to  the  Y  axis  of  a 
Hewlett-Packard  model  7001  X-Y  plotter.  The  Triad 
G-10  transformer  between  the  bridge  and  preamplifier  is 
used  with  junctions  having  a  resistance  of  200  0  or  less  at 
1000  Hz.  The  bridge  and  analog  recording  system  functions 
in  the  same  way  as  that  described  in  I  with  a  slightly 
better  signal-to-noise  ratio.  The  modulation  supply 
(Fig.  2)  is  essentially  that  described  in  I.  The  sweep 
circuit  (Fig.  2)  is  similar  to  that  in  1  but  has  been  changed 
so  that  to  allow  convenient  positive-,  negative-,  or  center- 
zero  sweeps  to  be  executed.  The  sweep  is  provided  by  a 
motor  drive  to  one  of  the  2  kfl  (10  turn)  potentiometers  in 
the  circuit.  Using  a  reference  frequency  of  1000  Hz  I  he 
recorder  plots  a  signal  proportional  to  dV/dl  as  a  function 
of  junction  bias  V.  With  a  reference  frequency  of  2000  Hz, 
cPV/dP  is  plotted  as  a  function  of  V.  The  system  is 
arranged  in  such  a  way  that  analog  and  digital  recordings 
are  made  simultaneously. 

The  analog-to-digital  conversion  is  accomplished  with  a 
Yidar  5203  data  acquisition  system;  however,  other  similar 
systems,  or  components  as  outlined  below  can  be  used. 
With  current  prices  dropping  one  can  duplicate  the  func¬ 
tions  shown  in  the  blocks  (Fig.  1)  at  somewhat  lower  cost. 
We  now  describe  the  digital  part  of  the  system.  Our  system 
has  a  10-channel  reed  relay  scanner  (Yidar  604).  The 
scanner  is  arranged  so  that  the  signal  from  the  lock-in 
amplifier,  our  dependent  variable  proportional  to  dV/dl 
or  cPV/dP,  is  fed  to  channels  2-5  and  7-10  in  parallel. 
The  junction  bias  voltage  is  our  independent  variable 
and  it  is  fed  through  a  buffer  amplifier  (Hewlett-Packard 
2470A)  to  channels  1  and  6  of  the  scanner.  Since  the 
scanner  operates  in  a  continuously  cyclic  mode  and  the 


Fig.  1.  The  bridge  circuit. 


Fig.  2.  The  sweep  circuit  and  the  modulation  supply. 


sweep  progresses  in  time,  we  obtain  a  value  of  bias  voltage 
V  by  averaging  channels  1  and  6.  The  corresponding  de¬ 
pendent  variable  (dV/dl  or  d-V/dP)  is  obtained  by 
averaging  the  output  of  channels  2-5.  The  next  pair  of 
points  is  obtained  by  forming  V  from  an  average  of 
channels  6  and  1  and  the  dependent  variable  by  averaging 
channels  7-10,  etc.  One  can,  of  course,  use  only  two 
channels  for  the  recording,  but  we  have  been  using  this 
averaging  procedure  because  transient  noise  from  other 
abs  in  the  building  is  sometimes  troublesome.  The  scanner 
signal  is  fed  to  an  integrating  digital  voltmeter  (\  idar 
521).  This  IDVM  is  usually  operated  with  a  166  msec  in¬ 
tegration  time  which  provides  a  resolution  of  one  part  in 
105.  The  coupler  (Vidar  653-03)  provides  the  interface 
between  the  IDVM  and  the  7-track  BCD  coded  tape,  it 
also  provides  the  requisite  command  signals  for  the  scanner 
and  voltmeter  to  read  the  data  sequentially.  Again  the 
recording  and  interfacing  can  be  done  with  other  compo¬ 
nents;  in  fact  9-track  tape  and  a  less  redundant  coding 
could  be  used. 

The  absolute  calibration  of  the  conductance  is  accom¬ 
plished  by  replacing  the  tunnel  junction  with  a  precision 
decade  resistance  and  recording  the  lock-in  output  for 
about  15  resistance  values  spanning  the  resistance  varia¬ 
tion  of  the  junction.  T  he  data  are  recorded  on  four  separate 
files  on  the  tape  to  allow  ease  of  computer  handling.  The 
first  file  contains  nominally  2000  pairs  of  points  pertaining 
to  the  bias  voltage  and  the  uncalibrated  conductance 
(lock-in  output).  The  next  file  contains  15  pairs  of  calibra¬ 
tion  points  for  calibrating  the  conductance.  I  he  third  file 
contains  about  4000  uncalibrated  </T  ,  dP  points  and  their 
respective  bias  voltages.  1  he  fourth  tile  contains  the  zeio 
level  for  the  calibration  of  the  second  derivative.  This  is 
obtained  bv  recording  the  lock-in  output  level  flock-in 
tuned  to  2000  Hz)  when  no  harmonic  (2000  Hz)  signal  is 
present. 
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Fig.  3.  The  flow  chart  for  the  complete  computer  program. 


AND  S  T  K  A  U S 

II.  CALIBRATION 

In  this  section  we  describe  the  basic  calibration  scheme 
used  to  obtain  values  of  a  and  da/ dV.T  he  actual  procedure 
and  flow  chart  for  programming  of  this  information  are 
given  at  the  end  of  this  section. 

A.  Conductance  (a)  Calibration 

Let  the  lock-in  output  at  the  fundamental  frequency  be 
denoted  as  y\.  Then  as  shown  in  I  the  junction  dynamic 
resistance  d  V / dl  is  related  to  Vi  in  accordance  with  Eq.  ( 1 ) - 
The  first  file  on  the  magnetic  tape  is  a  collection  of  2000 
points  (V,yi)  and  thus  contains  the  necessary  information 
to  obtain  a  plot  of  a  as  a  function  of  V  if  we  know  the 
calibration  of  Vi  in  terms  of  a  known  set  of  resistance 
values  R.  This  information  is  recorded  on  the  second  file 
and  because  the  bridge  response  is  linear  we  can  obtain  m 
and  b  for 

a  =  dI/dV  =  (myi+b)  \  (d) 

from  these  data,  thus  obtaining  a  calibrated  conductance 
plot  (o-  as  a  function  of  V ). 

B.  Second  Derivative  ( dv/dV )  Calibration 

Denoting  the  second  harmonic  output  of  the  lock-in  as 
y-i  we  have  from  1 

da/d  V = (PI /d  V2  =  rja3(V)y-i(V),  (4) 

where  y-iiV)  K  dPV  / dP. 

Because  of  small  amounts  of  second  harmonic  distortion 
present  in  the  modulation  supply  (0.01%),  we  cannot  treat 
rj  as  a  constant  (see  Appendix)  so  we  will  write  ??  =  i?(l  )• 
We  can  then  proceed  to  determine  tj(I')  from  <r(l  )  and 
y->(V)  using  the  integral  equation 

<r(V)-er(Va)=  j  V  (x)aM(x)yt(x)dx,  (•'>) 

J  V  a 

where  Va  is  any  particular  value  of  bias  voltage.  We  now 
write  Eq.  (5)  in  an  operator  equation  form, 

ku  =  f,  ((>) 

where 

f(V)=a(V)-a{Va).  (0 

If  Eq.  (6)  is  solvable,  then  the  solution  r?  will  minimize  the 
functional 

F(u)  =  \\Au-f\\2=  [  {Au-jfdx.  (S) 

J  Va 

The  replacement  of  Eq.  (6)  by  the  variational  problem  of 
minimizing  the  functional  F (u)  is  essentially  a  least  squares 
technique.10  We  can  expand 

n 

Un=H  (Ik'pk, 

k= t 


(9) 


- 

' 
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and  choose  the  expansion  coefficients  so  that  ||A«n  /I!2 

is  a  minimum.  Using  the  standard  notation  foi  the  scalar 
product  of  two  functions  f(x)  and  g(x)  as 

(/,£)  =  f  f(x)g(x)dx>  (10) 


we  obtain  a  system  of  n  linear  algebi  aic  equations 
(i  =  l,  2,  .  .  n), 


X:  ak(A\J/k,A\pi)  —  (A\pi,f).  (11) 

i=l 


Equations  (5)  and  (6)  can  now  be  put  in  this  format  choos¬ 
ing  the  set  {*pk}  in  polynomial  foim  \ pk  %  •  Then 


(Ax*-1, Ax*-1)  = 


•  V  |- 


-■>  v, 


X 


rz(v)yi(v)vl  ldv 


w  V  a 


dz  (12) 


and 

{Axi-\f)  = 


V  r  r  z 


J  Va 


az{v)yAv)v'  l(dv 


Va 


X[o(z)-o(Va)Jdz,  (13) 


where  i,k  —  1,  2  ...  n.  In  this  way  one  can  solve  Eq.  (11) 
and  obtain  a  polynomial  solution  for  77(F) 

r ;  (F)  =  6I1+O2F +aaF2+  •  •  • .  (14) 

One  finds  in  practice  that  a  quadratic  form  for  77  is  suffi¬ 
cient  to  determine  absolute  values  of  77  to  a  few  parts  in  103. 

We  now  outline  the  programming  scheme  to  implement 
the  above  method.  Figure  3  shows  a  flow  chart  for  the  com¬ 
puter  programming.  Program  parameters  and  options  are 
specified  in  block  1.  These  parameters  include  the  following 
information : 


(i)  number  of  files  to  be  read  on  input  tape; 

(ii)  range  of  the  independent  variable,  F(e.g.,  -300 

to  +300  mV) ; 

(iii)  interval  of  F  used  in  output  data  (e.g.,  a  point  at 
every  0.2  mV  voltage  change) ;  and 

(iv)  resistance  calibration  data,  i.e.,  values  of  resistance 
corresponding  to  the  Yi  values  used  in  calibration  (see 

Sec.  II. A). 

The  available  options  include  a  specification  of  which  and 
how  many  of  the  dependent  variables  are  to  be  calculated 
and  plotted.  The  plotting  is  done  by  generating  a  7-track 
BCD  tape  which  then  is  used  to  make  off-line  plots  on  a 
Calcomp  plotter.  Since  in  some  of  our  measurements  we 
are  interested  in  inelastic  processes  which  occur  in  tunnel 
junctions  and  these  are  even  functions  of  bias  voltage  F, 


it  is  convenient  to  separate  a  and  do/ dV  into  its  even  and 
odd  parts, 

oeiy)=\io{y)+o{-v)-]  (is) 

<r0(F)=40(F)-<r(-F)].  (16) 

Furthermore,  it  is  sometimes  convenient  to  reduce  these 
functions  to  dimensionless  quantities  independent  of  actual 
junction  resistance,  in  order  to  facilitate  intercomparison 
of  junctions.  We  achieve  this  by  arbitrarily  normalizing 
the  junction  conductance  by  dividing  it  by  its  value  at 
zero  bias,  i.e.,  y  —  o{V)/o(0).  Thus  in  block  1  we  choose 
from  the  following  dependent  variables:  0,  oe,  o0,  y,  ye,  y<>, 
do/dV,  doe/dV,  doo/dV,  dy./dV,  dyJdV,  dy0/dV. 

Block  2  in  Fig.  3  shows  the  conductance  data  input  (file 
1  on  data  tape).  Block  3  (file  2  on  data  tape)  includes  the 
input  data  for  calibration  and  calculates  m  and  b  in 
Eq.  (3).  The  deviations  from  linearity  are  also  calculated 
and  an  optional  plot  can  also  be  made.  Since  in  practice 
with  a  digital  computer  calculations  are  made  on  a  fixed 
grid,  i.e.,  equal  spacing  of  independent  variable,  the  2000 
or  so  points  (V ,Y \)  are  reduced  to  about  500  points  at 
fixed  voltage  intervals  (say  0.2  mV).  This  procedure  is 
represented  by  block  4  along  with  a  calculation  of  the 
corresponding  ( V,o )  points  using  Eq.  (3).  Blocks  5-7 


Fig.  4.  The  flow  chart  for 
the  integral  equation  [Eq. 
(5)3  subroutine. 
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Fig.  7.  The  odd  part  of  the  conductance  of  the 
junction  shown  in  Fig.  5. 
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FTg.  5.  The  junction  conductance,  normalized  by  the  conductance  at 
zero  bias,  of  a  Pb-insulator-Pb  junction  at  1.05  K  with  8  kOe  field. 


separate  the  conductance  into  odd  and  even  parts  and  cal¬ 
culate  the  required  junctions  and  produce  the  plotting  tape. 

We  now  turn  our  attention  to  the  calculation  of  absolute 
values  of  da/dV.  The  4000  or  so  pairs  of  points  (V,Y 2) 
(file  3)  and  the  zero  level  (file  4,  described  Sec.  I)  are 
the  inputs  in  block  8.  Block  9  compares  this  to  the  zero 
level  obtained  by  calculating  the  location  of  a  minimum 
in  the  a(V)  points  from  the  previous  calculations.  This 
provides  a  self-consistency  check  for  those  junctions  which 
have  conductance  minima.  In  order  to  carry  out  a  solution 
of  the  integral  equation  (5)  we  have  to  insure  that  we 
have  Y 2  values  at  the  same  voltage  points  for  which  we 
have  conductance  values.  This  overlap  and  interpolation 
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Fig.  6.  The  even  part  of  the  conductance  of  the 
junction  shown  in  Fig.  5. 


function  is  achieved  in  blocks  10  and  11.  The  solution  of 
the  integral  equation  (block  12)  is  shown  in  more  detail 
in  Fig.  4.  Following  the  input  of  data  the  matrix  elements 
are  calculated  in  block  2  with  w  =  3  (larger  values  of  n 
have  not  been  required  for  our  junctions).  A  solution  of 
this  system  of  three  linear  equations  in  three  unknowns 
provides  us  with  the  three  coefficients  ah  a2)  a3  of  Eq.  (14). 
In  practice  7?  is  calculated  separately  for  positive  and 
negative  voltages.  The  rest  of  the  calculations  calibrate 
da/dV  bv  Eq.  (4)  and  perform  a  check  by  calculating 

rv  da 

ct'(F)  —  o-(F.4)+  /  - d\  (1/) 

JradV 

and 

e(y)  =  la'(V)-a(y)ya(V)  (18) 

for  all  the  500  or  so  values  of  V  used,  t  is  always  within  a 
few  parts  in  105  (see  Appendix).  Returning  now  to  Fig.  3, 
block  13  calculates  odd  and  even  parts  of  da/dV 


As  a  double  check,  dae/dV  and  da„/d\'  are  integrated 
by  Eq.  (17)  and  compared  to  the  values  of  ae  and  <r„  ob¬ 
tained  directly  (block  5) ;  again  the  discrepancy  is  only  a 
few  parts  in  1()5.  The  rest  of  the  program  simply  outputs 
the  rest  of  the  data  and  plots. 


III.  OPERATION  AND  RESULTS 

The  bridge  is  balanced  at  some  value  of  bias  within  the 
sweep  range  in  such  a  way  that  the  full  dynamic  range 
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116  r 


08 


I  00 


0  92 


T  0  84 
> 

E  0  76 


o 

-  0.68 


0  60 


0  52 


044 


0  36 


0.28 


0  20 


1 


0  00  7.50  15.00  22  50  3000  3750  4500  52  50  6000  6750  7500 

ENERGY  (mV) 
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of  the  junction  shown  in  Fig.  5. 


of  the  lock-in  amplifier  (tuned  to  1000  Hz)  is  employed  in 
the  process  of  sweeping.  The  capacitance  Cd  is  adjusted 
to  compensate  for  the  junction  and  stray  capacitance. 
The  sweep  is  then  initiated  and  analog  and  digital  re¬ 
cording  is  carried  out  simultaneously.  Following  comple¬ 
tion  of  the  sweep,  an  end  of  file  mark  is  placed  on  the  tape. 
The  junction  is  now  replaced  by  a  decade  resistance,  Cd 
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Lio.  9.  The  even  part  of  the  derivative  of  the  conductance 
of  the  junction  shown  in  Fig.  5. 


is  removed,  and  calibration  points  are  recorded  on  the 
second  file.  Typical  results  for  a  Pb-insulator-Pb  junction 
are  shown  in  Figs.  5-7.  These  were  measured  at  1.05  K 
with  an  8  kOe  field  to  keep  the  Pb  from  becoming  super¬ 
conducting.  The  second  derivative  measurements  are 
carried  out  by  retuning  the  lock-in  to  2000  Hz  while  the 
bridge  is  still  fed  with  the  1000  Hz  signal.  Simultaneous 
analog  and  digital  recording  is  again  carried  out,  the 
digital  recording  being  placed  in  the  third  file  of  the  tape. 
Following  the  completion  of  the  sweep  the  lock-in  amplifier 
is  shorted  at  the  input  and  its  output  voltage  is  recorded 
in  the  fourth  file  to  provide  the  zero  level  reference.  Typical 
results  are  shown  in  Figs.  8-10. 

This  system  has  now  been  in  operation  for  about  two 
years  and  is  found  to  be  simple  and  convenient  to  operate. 
The  data  shown  in  the  figures  were  obtained  by  recording 
about  2000  first  derivative  points  and  4000  second  deriva¬ 
tive  points.  The  output  plots  for  y  and  dy/dV  (Figs.  5 
and  8)  contain  1000  points  while  the  remaining  odd  and 
even  plots  have  500  points.  The  computing  time  required 
was  2  min  on  an  IBM  360-67  computer,  using  Fortran  H. 
The  off-line  plotting  time  for  these  plots  requires  about 
15  min. 

We  wish  to' thank  Dr.  j.  S.  Rogers  for  many  useful 
discussions  during  the  course  of  this  work. 

APPENDIX:  EFFECT  OF  HARMONIC  DISTORTION 

Consider  the  simplified  bridge  in  Fig.  11.  Let  the  modu¬ 
lation  frequencv  be  co  and  allow  for  a  small  fraction,  .v,  of 
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second  harmonic  distortion.  We  will  neglect  higher  har¬ 
monics  in  this  analysis.  We  can  then  describe  the  current 
modulation  through  Rd  as 

I m  — 5 (cos>cot-\~ x  cos2cu/).  (Al) 

Note  x  is  a  very  small  fraction  in  our  bridge  0.0001. 
Let  the  tunnel  have  a  dynamic  conductance  aT=  (dl/dV) 
and  define 

r(V)=Rda(V)-l  (A2) 

T'(V)  =  Rdda/dV,  (A3) 

then  at  balance  we  have  T  =  0. 

For  simplicity  we  denote  the  modulation  voltage  ampli¬ 
tude  across  Rd  as 

td  =  Rd8.  (A4) 


Then  the  lock-in 
will  be 


Fu.cc 


output  at  the  fundamental 


r 


i'd 


r+i 


r  - 

r(r+i)2_' 


frequency 


(A5) 


The  second  harmonic  output  will  be 


'M2  r' 

Fjju  cc  l  ~  ) 


2/  (r+i): 


4x  (r+i)2- 
*d  T'  J 


(A6) 


We  note  that  the  fractional  error  in  the  first  derivative 
increases  linearly  with  increasing  modulation  level  while 
the  second  derivative  fractional  error  decreases  with  in¬ 
creasing  modulation. 

If  the  bridge  is  balanced  and  a  linear  element  is  put  in 
the  place  of  the  junction  (da/dV  =  0),  then  the  harmonic 
distortion  has  no  effect  since  it  cancels  on  both  sides.  In  the 


Fig.  12.  The  junction  conductance,  normalized  by  the  conductance 
at  zero  bias,  of  a  Co-CoO-In  junction  at  4.2  K. 

case  of  a  nonlinear  element  (da/dV^O),  however,  the 
effect  is  to  unbalance  the  bridge  thereby  preventing  can¬ 
cellation  for  the  small  distortion  signal.  This  is  the  reason 
that  r]  is  not  a  constant  and  has  the  form  given  in  Eq.  (14). 
For  the  most  part  calibrating  a  by  substituting  a  resistance 
decade  along  with  the  computations  outlined  in  Sec.  11 
provides  a  good  calibration  with  e  [Eq.  (18)]  kept  to  a 
few  parts  in  105.  Very  nonlinear  junctions  such  as  lig.  12 
will  give  e~  10~3.  But  for  these  junctions  we  carry  out  the 
recording  piecewise  so  that  Rd  is  readjusted  various  times 
in  the  sweep  range  to  rebalance  the  bridge;  in  this  way  t 
can  again  be  reduced  to  a  few  parts  in  105. 
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The  electron  bottleneck  in  normal  metal— insulator— metal  tunnel  junctions 
due  to  nonequilibrium  effects  has  been  investigated  in  the  presence  of  a 
magnetic  field.  The  observed  results  are  due  to  a  new  geometrical  size  effect 
and  are  explained  in  an  extension  of  our  nonequilibrium  theory. 


IT  HAS  recently  been  shown1  that  small  conductance 
minima  (a  =  d//d  V)  below  the  phonon  structure 
observed  in  normal  metal— insulator— metal  junctions 
are  due  to  the  finite  relaxation  times  of  the  tunneling 
electrons  in  the  metals  comprising  the  junction.  States 
which  would  otherwise  be  available  for  tunneling  are 
blocked  creating  an  electron  bottleneck.  This  non¬ 
equilibrium  effect  has  been  explained  in  detail  in  a 
transport  model. 

In  the  previous  work,1  superconductivity  in  the 
electrodes  of  the  junction  was  quenched  using  a 
magnetic  field.  In  this  letter  we  report  on  a  detailed 
study  of  the  magnetic  field  dependence  of  the  small 
Zero-Bias-Anomalies  (ZBA)  observed  in  such  tunnel 
junctions.  An  explanation  is  given  for  the  observed 
behavior  as  a  size  effect  due  to  the  finite  geometry 
of  the  junction. 

For  a  systematic  study  of  magnetic  field  effects 
in  tunnel  junctions  it  is  useful  to  think  in  terms  of 
two  possible  field  orientations:  the  longitudinal  field 
which  is  parallel  to  the  tunneling  electron  momentum 
thus  providing  no  Lorentz  force  on  the  electron,  and 
the  transverse  field  which  is  perpendicular  to  the 
tunneling  electron  motion  giving  rise  to  a  Lorentz 
force.  Measurements  on  a  number  of  junctions  using 
different  electrode  metals  (Al,  Au,  Mg,  Pb,  Bi)  have 
shown  no  magnetic  field  dependence  of  the  ZBA  in 
longitudinal  fields  up  to  70  KG  (the  highest  field 

*  This  work  was  supported  in  part  by  grants  from  the 
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available  to  us).  This  suggests  that  in  the  low  energy 
region  where  the  ZBA  is  observed  no  appreciable  changes 
in  the  Fermi  surface  occur  in  our  junctions  that  could 
affect  the  tunneling  mechanism.2  However,  in  a 
transverse  field  a  dramatic  field  dependence  is  observed 
as  is  shown  in  Fig.  1(b)  which  gives  a  family  of  o(V ) 
vs  V  curves  as  a  function  of  magnetic  induction  B. 

For  comparison  Fig.  1(a)  shows  the  temperature  de¬ 
pendence  of  the  same  junction.3 


FIG.  1.  Normalized  even  part  of  the  conductance  oe 
(F)/oe  (V  —  5  mV)  vs  V  for  an  Al— Al  junction  (85  1  A) 
for  (a)  different  temperatures  T  and  (b)  various 
magnetic  inductions  B. 


In  our  previous  work1  we  have  shown  that  the 
electron  bottleneck  arises  from  the  competition  of  the 
tunneling  rate  D/t0  with  the  elastic  scattering  frequency 
t,'1  (e.g.  impurity  and  diffuse  boundary  scattering)  and 
the  electron-phonon  relaxation  rate  (co,  T ),  the 
latter  giving  rise  to  the  temperature  dependence  of  the 
ZBA. 
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To  understand  the  magnetic  field  dependence  we 
need  not  consider  the  electron-phonon  interaction 
because  as  mentioned  above  Fermi-surface  effects  are 
of  no  consequence  here.  Furthermore,  because  the 
tunneling  barrier  is  less  than  50A  thick  the  tunneling 
rate  will  be  unaffected  by  the  Lorentz  force.  However, 
the  Lorentz  force  induced  electron  drift  in  the  elec¬ 
trodes  provides  an  additional  mechanism  to  take 
electrons  out  of  the  forward  direction  thus  expanding 
the  bottleneck  by  providing  more  tunneling  states 
in  the  electrodes.  In  the  absence  of  a  magnetic  field 
specular  scattering  of  electrons  by  the  electrode 
boundaries  will  not  be  effective  in  relieving  the 
bottleneck  because  such  scattering  does  not  remove 
electrons  from  the  forward  direction.  In  a  transverse 
field,  however,  the  specularly  reflected  electrons  will 
drift  along  the  electrodes.  This  provides  an  additional 
mechanism  for  the  removal  of  electrons  from  the 
forward  direction.  In  our  transport  model  we  can 
account  for  this  drift  by  a  relaxation  time 


rd{B)  = 


a 

‘'drift 


a  l 
vp  R 


,2  V/21  _1 

*7  - 


where  a  is  the  effective  width  of  the  junction  in  the 
direction  of  the  drift,  vF  is  the  Fermi  velocity  of  the 
electron,  /  is  the  electrode  thickness  and  R  =  mvpc/eB. 
For  B  <  60  KG  R>  20, 000 A  which  is  much  larger 
than  the  electrode  thickness  used  here  (/  <1  1000 A). 
Hence  we  can  simplify,  rd(B )  ~  2 mc/e'd/l' B~x  and 
the  total  relaxation  time  in  one  electrode  is  given  by 
r'1  =  tJ1  +  rdx  C B )  +  t;x(co,  T).4  To  simplify  the 
discussion  let  us  assume  that  the  bottleneck  is 
entirely  due  to  one  electrode  only.  (This  can  be 
achieved  e.g.  by  using  a  second  electrode  with  a  strong 
electron— phonon  interaction.)  We  then  find  for  T  = 
OK  and  assuming  that  the  electron-phonon  coupling 
is  linear  in  co  (a2F  =  a2  to)  in  the  low  energy  region, 


o(V)  «  D 


D/t0 


l/r,+  '-B  +  l  a3(eV)2 
2  me  a  h 


(1) 


This  shows  that  the  effect  of  B  is  to  reduce  the  ZBA 
by  removing  electrons  from  the  forward  direction. 
The  rate  at  which  this  reproval  takes  place  for  a  given 
B  depends  entirely  on  the  electrode  geometry  via  the 
ratio  l/a  of  the  electrode  thickness  /  to  the  effective 
width  of  the  junction  a.  This  size  effect  has  been 
checked  experimentally  by  varying  the  film  thickness 


l  and  by  looking  at  crossed  film  junctions  in  which 
the  width  of  base  and  top  films  were  different,  thus 
giving  rise  to  a  rectangular  tunneling  area  of 
dimension  ax  •  ay  with  ax  =£  ay.  It  was  found  that  a 
predictable  change  occured  according  to  equation  (1) 
when  the  field  was  changed  from  the  x -direction  into 
they-direction. 
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FIG.  2.  (a)  Peak  heights  (open  circles)  and  peak  positions 
(full  dots)  in  the  normalized  derivatives  of  the  even 
part  of  the  conductance  as  a  function  of  temperature  for 
an  Al-Al  junction  (851  A);  triangles;  data  for  the 
691  A  junction,  (b)  Peak  heights  and  peak  positions 
for  the  85 1  A  junction  as  a  function  of  magnetic 
induction  B.  Solid  curves:  theoretical  results. 


We  have  analysed  quantitatively  the  data  for  an 
Al-Al  junction.  Fig.  2(a)  shows  the  peak  heights  and 
peak  positions  in  the  normalized  even  derivative  of  o(V), 
i.e.  l/o(Q)'doe(V)/dV  as  a  function  pf  temperature. 

To  show  the  reproducibility  of  our  data,  we  have 
included  in  this  graph  data  taken  on  a  different 
junction  a  year  ago.1  (Note:  A  direct  comparison 
of  the  two  junctions  can  only  be  made  for  the  peak 
positions  which  depend  on  t,  and  rep  alone.  Without 
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renormalization  a  comparison  of  peak  heights  is  not 
meaningful  because  the  peak  heights  are  proportional 
to  the  (geometry  dependent)  tunneling  rate  D/t0  which 
differed  by  about  a  factor  of  5  for  the  two  junctions 
due  to  slightly  different  resistances  and  much  thinner 
A1  films  in  the  junction  851.)  The  elastic  scattering 
time  obtained  from  the  theoretical  curves  is  r,  =  2.5. 10" 10 
sec  assuming  a2  =  0.003/meV  for  A!.1 

To  make  a  similar  analysis  for  the  field  dependence 
we  can  use  equation  ( 1 )  for  T  =  OK  to  obtain 

1  da(Fmax,  5) 

o(0,  5)  dV 

1  do(f  max,  0)  C'T,  m-3 

o(0, 0)  dF  2 me  a  2  (2) 

Fma x(B)  =  ^max(9)  ( I  +  ^  ^ 

for  the  variation  of  the  peak  height  and  peak  position 
in  da/dF  as  a  function  of  B.  It  should  be  noted  that 
equations  (2)  are,  of  course,  only  rigorously  valid  for 


T  =  OK.  But  at  low  enough  constant  temperature  the 
5-dependence  should  be  fairly  accurately  described. 

The  theoretical  curve  in  Fig.  [2(b)]  was  obtained  by 
using  the  value  r,-  =  2.5  X  10"10sec  from  the  temperature 
dependence  and  the  geometrical  data  of  our  junction, 
namely  l  =  680  A  and  a  =  0.3  mm.  This  leads  to  a 
coefficient  5  X  10”3/kG  in  front  of  B  in  equations  (2), 
which  reproduces  the  data  surprisingly  well.  It  should 
be  noted  once  more  that  the  theoretical  fit  of  the 
5-dependence  is  parameter-free.  The  quality  of 
agreement  might  be  viewed  as  fortuitous  but  certainly 
an  order-of-magnitude  agreement  should  always  be 
possible. 

In  conclusion,  the  magnetic  field  dependence  of 
the  small  ZBA  in  the  transverse  case  has  been  explained 
in  terms  of  a  size  effect  due  to  the  finite  geometry 
of  the  tunneling  junction.6  Moreover,  the  absence  of 
any  field  effects  on  the  ZBA  in  the  longitudinal  case 
has  been  understood  and  leads  to  a  clear  separation 
of  the  non-equilibrium  effects  at  low  bias  from 

<y 

possible  Fermi-surface  effects. 
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5.  JAKLEVIC  R.C.,  LAMBE  J„  MIKKOR  M.  and  VASSELL  W.C.,  Phys.  Rev.  Lett.  26,  88  (1971);  in  this 
reference  the  ‘quantum  box’  size  effect  is  investigated  which  is  quite  different  from  the  size  effect  considered 
here  which  is  connected  with  the  electron  trajectories. 


Der  Einfluss  eines  Magnetfeldes  auf  gewisse  Transportphenomene  beim 
Elektrontunneleffekt  in  Metall-Isolator-Metall  Tunneldioden  wird 
untersucht.  Die  Resultate  werden  als  ein  neuer  geometrische  Effekt  gedeutet 
und  in  einer  Erweiterung  unserer  Nichtgleichgewichtstheorie  quantitativ 
erklart. 
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